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1 Introduction 



Quantum computers p][2l[3lH] are explicitly quantum mechanical systems that 
use phenomena such as superposition and entanglement to perform computa- 
tional tasks more efficiently than any classical computer [5j. More generally, 
quantum information technologies can perform information processing tasks, 
such as unconditionally secure communication, which cannot be achieved by 
classical systems. Unsurprisingly quantum computation has generated an enor- 
mous amount of interest, reflecting not just its potential technological impor- 
tance, but also the intellectual importance of the challenge it provides to 
previous formulations of computational complexity theory [6], built on the 
Church-Turing thesis, which asserts that all physical models of computation 
are essentially equivalent to one another in computational power. 

This intrinsic interest is partly tempered by the great practical difficulty in 
building large scale devices capable of implementing technologically impor- 
tant computations, but it has at least proved fairly simple to build small 
demonstration devices, and NMR has played a leading role in this field for 
many years. The ffist ideas on how to build quantum computers with NMR 
[7f5|5piipT| were swiftly followed by the first implementations of quantum 
algorithms p^|T3] , and it is only recently that other techniques have achieved 
similar sophistication in the manipulation of quantum states: although an im- 
plementation of the fundamental controlled-NOT gate with trapped ions was 
demonstrated in 1995 [14j it was not until 2003 that a complete implementa- 
tion of Deutsch's algorithm was achieved [I5j. It must, however, be remem- 
bered that the great difficulty in preparing NMR systems in pure spin states, 
reflecting the tiny energy gap between nuclear spin states, has given rise to 
grave concerns about the direct relevance of NMR techniques to attempts to 
build large scale devices [ISIIZ], and has even led to questioning of whether 
many NMR quantum computations can be considered true quantum compu- 
tations at all [18] . 

1.1 Structure and scope 

In my previous introduction [19] and reviews [20|21|22] I described the funda- 
mental concepts of quantum computation, and how these are related to more 
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familiar concepts in conventional NMR spectroscopy. For the most part I will 
only briefly repeat these here, and will largely assume familiarity with these 
basic concepts. A large number of excellent textbooks, e.g. |2H|24II25|I2B] . have 
now been published, providing many different approaches to the necessary 
background with almost any desired style and level of sophistication, and a 
number of reviews, e.g. |27ll28|29|30|31|32j . have described both fundamen- 
tal concepts and recent developments in particular areas of NMR quantum 
computing. 

The general requirements for implementing quantum information processing 
on a physical system are usually considered in terms of the seven DiVincenzo 
criteria [33fM] , where the first five criteria are required for quantum computa- 
tion and the final two permit quantum computers to be connected by quantum 
information links, in effect building a quantum internet. The first five criteria 
are: 

(1) a scalable physical system with well characterized qubits; 

(2) the ability to initialize the state of the qubits to a simple fiducial state 
such as |000 . . . ); 

(3) long relevant decoherence times, much longer than the gate operation 
time; 

(4) a "universal" set of quantum gate; 

(5) a qubit-specific measurement capability. 

As is traditional I will structure the first part of my discussion around these 
five criteria, but will take them in a different, and more helpful, order. I will 
not consider the two additional criteria, the interconversion of static qubits 
(suitable for information processing) and flying qubits (suitable for information 
transmission), and the transmission of these flying qubits from place to place. 

I will only describe in detail the simplest method for implementing quantum 
information processing with NMR, that is using spin-| nuclei in molecules in 
solution, although I will briefly comment on alternative approaches including 
solid-state NMR, the use of liquid crystal solvents, and the use of quadrupo- 
lar nuclei, as well as on approaches involving electron spins, either alone or 
combined with nuclear spins. As the discussion is structured around the Di- 
Vincenzo criteria, these topics will inevitably be scattered around the text. I 
will not discuss exotic long range proposals, such as the famous Kane proposal 
[35] for a large scale quantum computer based on ENDOR manipulation of '^^P 
spins in doped silicon, but will largely confine myself to methods which can 
be implemented with relatively conventional NMR and ESR spectrometers, 
and for which at least elementary experimental demonstrations of quantum 
information processing have in fact been performed. These sections are fol- 
lowed by an introduction to some quantum algorithms and communication 
protocols, highlighting some experimental achievements. The available litera- 
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ture is now too large to include every paper, or even every topic that has been 
investigated, and my selection is inevitably personal and incomplete. 



2 Qubits 

A quantum bit, or qubit [36], is simply a two level quantum system, where the 
two levels are labelled as |0) and effectively identifying them with the two 
logic states and 1. This choice is usually called the computational basis, and 
is the implicit basis used here; the state |0) is usually chosen to have lower 
energy than |1) but this is not essential. 

2.1 Superpositions and entanglement 

As a qubit is a quantum system it is not confined to these two basis states, 
but can exist in a general superposition state of the form 



where a and /3 are arbitrary complex numbers subject to the constraint that 
lap + = 1. However, this traditional description of a superposition state 
of a qubit is in some sense excessive, in that there is an ambiguity in the global 
phase of the state [37], so that the set of states 



(where 7 is the global phase, taking real values) are completely indistinguish- 
able from one another. Equivalently, the global phase can be incorporated into 
the value of one of the coefficients, usually taken as the first coefficient a, so 
that this coefficient can be assumed to be real. Such global phases must be 
sharply distinguished from the relative phase of two coefficients in a superpo- 
sition, which is a measurable quantity. 

The ability of quantum computers to process data in superposition states lies 
at the heart of their power. One particularly important pair of superposition 
states is |-|-) and |— ), defined by 



which play a key role in many quantum algorithms. 

A single qubit is not a particularly interesting object, but even a two qubit 
system can show interesting computational behaviour. A pair of qubits has 



a|0) + 





(2) 
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four basis states, and can be found in some general superposition of all four, 



a|00) + /3|01)+7|10) + 5|ll) (4) 

with the constraint |ap + + I7P + |5p = 1. These general states can 
be divide into two major groups: product states, which can be written as the 
direct product of two terms of the form Eq. [H such as 

. 10)1^). (5) 

and entangled states, which cannot be decomposed in this way, such as the 

St cXt cs 

|oo)±|ii) |oi)±|io) ... 

10.) = ^^ 1^.) = ^^ (6) 

which are the four maximally entangled Bell-states. Entanglement is an impor- 
tant, but surprisingly slippery, concept which will be returned to on a number 
of occasions. The slipperiness arises from the fact that although it is relatively 
easy to demonstrate that a state is a product state, by writing it down in 
explicit product form, it is much more difficult to prove that a state cannot 
possibly be written in a product form, and so apparently entangled states may 
in fact turn out to be product states. 



2.2 Pure states and mixed states 



A qubit in a well defined quantum state (which can be a computational 
basis state or a superposition state) is most conveniently described in the ket 
notation used above, but such a state, called a pure state, can also be described 
using the density matrix description 



I a\ / \ / aa* a/3* 



(7) 



\/3a* I3l3*^ 

where the constraint on a and /3 now appears as the property that 

tr(p) = |a|2 + |/3|2 = l. (8) 
For a pure state the density matrix p has the property that 

p' = mm){i^\ = m{i^\ = p (9) 

so that p is said to be idempotent, and tr(p^) = 1. A qubit can also be in a 
more general mixed state, which is a statistical mixture of pure states 

Ep#*)(V^*I (10) 
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where the Pi > are probabihties, so that J2iPi = 1- Density matrices for 
mixed states are not idempotent. Mixed states can be classified according to 
their purity, which can be defined as tr(p^) < 1, although other definitions are 
in use. Unitary transformations of a state leave the purity unchanged, while 
non-unitary operations normally act to to reduce the purity; increasing the 
purity of a state is a particularly tricky operation. 

The ambiguous global phase term is not present in the density matrix descrip- 
tion of a quantum state: for the state considered above, Eq. [21 the correspond- 
ing density matrix is 

e^^|^)(^|e-'^ = |^)(^| (11) 

as the contributions from the global phase terms in the ket and bra cancel out. 
Thus the density matrix contains only the physically observable information 
about a state. 

For a pure state, the ket and density matrix approaches give two different 
ways of relating qubit states to the Bloch sphere. Firstly, considering the two 
constraints on the state of a qubit, Eq. [H it can be written as 

cos(^/2)|0) + sin(^/2)e^'^|l) (12) 

where < < vr, which determined the relative amplitude of the two basis 
states, is a co-latitude, and < < 2tt, the relative phase, is an azimuth 
angle. These two angles define a unique point on a sphere of radius 1, the 
Bloch sphere. Alternatively, the density matrix description of a qubit, Eq. [7] 
can be written in the Pauli basis (which extends the three traditional Pauli 
matrices into a matrix basis by including (Tq, the two by two identity matrix) 
as 

I {SoO-Q + SxCTx + SyO-y + SzCTz) ■ (13) 

As both the density matrix and the Pauli basis are manifestly Hermitian, the 
four coefficients Sq, s^, Sy and are necessarily real. The constraint jap + 
I^P = 1 sets So = 1, leaving the remaining three coefficients which define a 
vector of length 1. A similar approach can be used for mixed states, but in 
this case the state is described by a point inside the Bloch sphere. This is 
most easily seen by noting that the Bloch vector for a mixed state will be 
a weighted vector sum of the contributing Bloch vectors, and so must have 
length less than one. 

In the discussions above I have used an explicit description of a density matrix, 
Eq. [ini composing it as a particular mixture of pure states. It is tempting to 
conclude that a given decomposition has some particular meaning, but in 
fact mixed state density matrices do not have unique decompositions, and no 
decomposition is more correct than any other. One famous example is that an 
equally weighted mixture of |-|-) and |— ) states is indistinguishable from an 
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equally weighted mixture of |0) and |1) states: 
ll+)(+l + l|-)(-| = 




(14) 



||0)(0| + i|l)(l| 



As the density matrix description contains all physically accessible properties 
of a state, these two mixtures are essentially equivalent. 

The belief that one particular decomposition of a mixed state has particular 
validity is sometimes called the preferred ensemble fallacy or partition ensem- 
ble fallacy [38], and can lead to serious confusion. This does not, of course, 
mean that it is not permissible to use a particular decomposition; indeed, it 
is frequently very useful to do so. However, one must be extremely cautious 
when making arguments that depend explicitly on any particular decomposi- 
tion. For example, it is tempting to argue that a mixture of entangled states 
must itself be entangled, but it is easily shown that the apparently entangled 
mixture 

+ (15) 
can also be decomposed as the mixture 

i|00)(00| + i|ll)(ll|, (16) 

which is clearly a mixture of product states. Such states are not entangled in 
any useful sense; this point will be explored in more detail when I consider 
the initialisation of quantum states. 



2.3 Spin-\ nuclei in liquid samples 



There are many possible physical implementations of a qubit, frequently based 
on a two-level subsystem of a larger quantum system, but a particularly nat- 
ural implementation is provided by a spin-| atomic nucleus, such as ^H, and 
this is the approach used in the great majority of NMR implementations. A 
two-qubit quantum computer can be built from two spin-| nuclei, and so on. 
As will be discussed below it is necessary that the two nuclei are distinct, so 
that the two qubits can be separately addressed, and there must be some sort 
of spin-spin interaction, so that two-qubit logic gates can be constructed. This 
is easily achieved by using two inequivalent nuclei in a molecule. The great 
majority of experiments have been performed using conventional liquid state 
samples, where rapid molecular tumbling greatly simplifies the spin Hamilto- 
nian [39] , with the isotropic part of the scalar J-coupling providing a suitable 
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form of coupling. Such a sample contains many copies of each molecule, but 
as intermolecular interactions are cancelled out by tumbling the result is an 
ensemble of identical independent molecules, which can for most practical 
purposes be treated as a single molecule in a mixed spin state. 

In systems with two or more qubits an important practical distinction can be 
made between systems where all the spins are different nuclear species (a fully 
heteronuclear spin system) and those with two or more nuclei of the same 
type (a partially homonuclear spin system). As there are a limited number of 
different spin-i nuclei, among which only six (^H, ^^C, ^^N, ^^F, ^^Si and ^^P) 
have been used in quantum computing experiments it is clear than only very 
small quantum computers can be fully heteronuclear, but the relative ease of 
working with such systems makes them popular for implementing simple tasks. 
Other spin- 1 nuclei, such as ^^^Sn and ^^^Sn could in principle be used, but 
have not yet found use, presumably reflecting their low gyromagnetic ratios, 
low natural abundances, and the relative difficulty in obtaining isotopically 
labelled compounds. Similarly the difficulty of obtaining NMR probes capable 
of addressing more than five nuclei simultaneously (including the lock), 
and the ready availability of sensitive probes optimised for particular nuclear 
combinations (such as ^H, ^'^C and ^^N, used in NMR studies of proteins) has 
to some extent determined the choice of spin systems used so far. 

A (very incomplete) list of systems of spin-| nuclei which have been used 
to implement NMR quantum computing is given in table [TJ Spins used to 
study single isolated qubits are not listed; nor are systems involving high-spin 
nuclei. A variety of larger spin systems have been explored in histidine [IQ] 
which contains 14 spin-| nuclei. 
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Table 1: A partial list of spin systems of spin-i nuclei 
used to implement quantum computing with NMR; only 
very selective references are given for each spin system. 
Isotopic labelhng is not indicated, and in some cases dif- 
ferent sized spin systems are created either by selective 
isotopic labelling of a compound, or by using decoupling 
to effectively remove a spin from the spin system. 



n 


spins 


molecule 


references 


2 


HH 


2,3-dibromothiophene 


[TiHnnEiis] 






cytosine 


[T21I53E3] 






uracil 


m 






5-nitrofuraldehyde 


piiTlliS] 






coumarin 


m 






Ir(H)Cl(H)(CO)(PPh3)2 


m 






Ru(H)2(CO)2(dppe) 


[50p] 






Ru(H)2(CO)2(dpae) 


[52] 






5-bromothiophene-2-carbaldehyde 


[53] 




HC 


chloroform 


[T3ll5ll55ll56llFfll5TO 






formate 


[601I6T1I62] 






benzene 


[63] 






dimethylformamide 


m 




HF 


5-fluorouracil 


[65] 




HP 


phosphonic acid 


[661167] 




PP 


1 ,4-diphosphafulvene 


[6H] 


3 


HHH 


2,3-dibromopropionic acid 


[69170] 






chlorostyrene 


[7T1I72] 




HHF 


4-fluoro- 7-nitro-benzof uran 


[^61173] 




HHP 


£'-(2-chloroethenyl)phosphonic acid 


m 




HCC 


trichlororethene 


[7511761177] 






tris(trimethylsilyl)silane- acetylene 


m 




HCF 


dibromofluoromethane 


m 



10 







diethyl- fluoromalonate 


m 




ccc 


alanine 


[771l8T|82ll83ll8i1l85l^^ 




FFF 


bromotrifluoroethene 








2,3,4-trifluoroaniline 


m 






iodotrifluoroethene 


[97198] 




PPP 


3 ,4-dihydro-l ,3 ,4-triphospliacyclopenta[ a] indene 


[68] 


4 


HHHH 


l-cliloro-2-nitrobenzene 






HHFF 


2 , 3- difluoro- 6-nit rophenol 


[1811731199] 




HCCC 


alanine 


m 




HNCC 


glycine 


[100] 




CCCC 


crotonic acid 


[T0T1IT02] 




FFFF 


3-cliloro-2 ,4, 5 , 6-tetrafluoropyridine 


[103] 


5 


HHFFF 


2,4,5-trifluorobenzonitrile 


|l()3fl()4] 




HNCCF 


glycine fluoride 


[105] 




HFFFF 


2,3,4,6-tetrafluoroaniline 


|l()f>fl()7fl()8] 




FFFFF 


pentafluorobutadienyl complex 


[109] 


6 


HHHHHH 


inosine 


m 


7 


HHHCCCC 


crotonic acid 


|lllfll2] 




CCFFFFF 


pentafluorobutadienyl complex 


m 
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One complexity which arises in larger spin systems is that they frequently 
include groups of two or more equivalent nuclei, such as the three nuclei in 
a methyl group, and these cannot be treated as independent qubits, but must 
all be considered together. However, it is possible to decompose a methyl group 
as the combination of a fictitious spin-| particle and a fictitious spin-| particle, 
and in some experiments it is possible to use only the spin-| component of a 
CH3 group [11111112] or the spin-1 component of a CH2 group 00]. Equivalent 
nuclei can also be used directly in situations where only partial control of the 
spin system is needed, such as entanglement assisted magnetic field sensing 



2.4 Product operators and deviation density matrices 



Liquid state NMR experiments are usually described using product operator 
notation [116] , which is almost equivalent to the Pauli basis, differing only in 
normalisation. As NMR spin states are almost always highly mixed it is com- 
mon practice to be quite casual about normalisation, and simply concentrate 
on the relative size of terms in the density matrix which can be measured. This 
approach can lead to serious errors when describing quantum computing, and 
it is necessary to proceed with care. 

We can describe single qubit density matrices using product operators; for 
example the pure ground state of a single qubit is described by 



|0)(0|= f^°l = f^°l + f^ ^]=IE + I^ (17) 



while |1)(1| = ^E — Iz- As NMR observables are traceless the component is 
not directly detectable, and so we can say that |0)(0| is equivalent to J^. This 
description, in which multiples of the density matrix are ignored, is sometimes 
called the deviation density matrix |TT] description of a state. Note that devi- 
ation density matrices are not proper density matrices, in that they have trace 
and negative eigenvalues, and so cannot actually describe the state of a real 
physical system; they are simply a useful computational fiction. Furthermore, 
the thermal state of a single spin is in fact the mixed state ^E + 5Iz, where S is 
the polarisation of the spin. Although the absolute value of 6 is not normally 
directly measurable, it will nevertheless ultimately determine the signal size, 
and so cannot be entirely ignored. 

For a single spin the thermal state is broadly equivalent to the desired initial 
state |0)(0|, but this simple equivalence does not remain true for larger spin 
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systems. For example, with two spins the desired initial state is 

/^l 0^ 



\^o oy 



|00)(00| 



-(- 



while the thermal state for a homonuclear spin system is \E + 5{Iz + Sz)- 
These points are explored in detail in Section H] where methods of initialising 
spin systems are considered. 



2.5 Dipolar couplings 



NMR studies in the solid state |117|118] are considerably more complicated 
than those in the liquid state as the direct dipolar coupling term in the spin 
Hamiltonian is not averaged by rapid molecular tumbling. This has the poten- 
tial advantage that direct dipolar couplings are normally considerably larger 
than the isotropic part of the scalar coupling, but this on turn means that 
dipolar couplings are unlikely to be in the weak coupling limit, complicat- 
ing the Hamiltonian further. More seriously, intermolecular dipolar couplings 
mean that it is not usually possible to treat the molecules as an ensemble 
of identical independent copies. Further difficulties arise from the orientation 
dependence of chemical shifts and the presence of multiple molecules with 
different orientations in the unit cells of many crystals. Nevertheless, the con- 
siderable advantages that solid state systems might offer [28j, most notably 
the abihty to cool nuclear spins close to their ground states, have encouraged 
initial attempts at this problem. 

The first implementation of quantum, computing by solid state NMR |119j 
used a single crystal of glycine, containing 10% of fully ^'^C and ^^N labelled 
glycine diluted with natural abundance glycine to minimise the effects of in- 
termolecular couplings. The sample orientation was set so that one molecule 
in the unit cell gave sharp lines well resolved from the signals from other orien- 
tations, and decoupling was applied throughout all experiments, effectively 
removing the nuclei from the spin system with the exception of an initial 
cross-polarization phase to increase the polarization of the ^^C and ^^N nu- 
clei. The heteronuclear couplings between ^^C and ^^N are weak compared to 
the Zeeman energy scale, and the homonuclear coupling between the two ^^C 
nuclei (conventionally called Co and Co) can also be treated (to a reasonable 
approximation) as weak. Thus the final spin system is very similar to that 
used in liquid state experiments, and similar techniques can be applied [119j . 
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This early work overcame the complexities inherent in solid state NMR by 
carefully choosing a spin system where they could be largely ignored rather 
than by tackling them head on. A later more ambitious approach |12Ufl21fl22] 
used three ^^C nuclei in fully ^^C labelled malonic acid. Once again an ori- 
ented single crystal was used containing 3.2% labelled malonic acid diluted 
with natural abundance malonic acid. In this spin system, however, the dipo- 
lar couplings cannot be taken as weak, and it is necessary to consider the full 
dipolar Hamiltonian in designing logic gates (the couplings are weak enough 
that it is still possible to identify the nuclear eigenstates with the computa- 
tional basis states in a reasonably straightforward manner). This problem was 
tackled by using strongly modulated composite pulses [89] and GRAPE (gradi- 
ent ascent pulse engineering) pulses (123] which will be described in Sections 
ESI and [321 

The methods described above all use oriented single crystals, as couplings 
which are truncated by the Zeeman interaction depend strongly on the rela- 
tive orientation of the coupling tensor (aligned with the crystal axes) and the 
magnetic field. The approach commonly used in conventional solid state NMR 
of spinning out these interactions and then selectively reintroducing them by 
recoupling |118] will not work well in NMR quantum computing as the orienta- 
tion dependence of the reintroduced couplings, visible as spinning sidebands, 
will prevent the accurate implementation of quantum gates. A radically dif- 
ferent approach, representing quantum information using Floquet states in 
spinning samples has been proposed [124j . but has not been further explored. 

NMR studies of molecules dissolved in liquid crystal solvents |125|126] pro- 
vide a convenient half-way house between conventional solutions and solid 
state samples: the partial local ordering induced by the solvent results in in- 
tramolecular dipolar couplings being scaled down, but not averaged to zero, 
while intermolecular couplings are effectively averaged, as in true liquids, and 
can be ignored. Thus the system can again be treated as an ensemble of iden- 
tical independent molecules. Early attempts to implement NMR quantum 
computation with liquid crystals used systems where the internal couplings 
could still be treated as weak [56|57fl03] . and the experiments are very simi- 
lar to conventional liquid state experiments, but later experiments have used 
systems with strongly coupled spins |127fl28] . 

A more extreme approach is to use clusters of equivalent spins; the dipolar 
coupling between such spins cannot normally be seen in liquids, but is revealed 
in liquid crystal solutions |129] . Early experiments were performed on the 
three nuclei and single ^^C nucleus in a methyl group |13U|131j . but these 
have now been extended to much larger spin systems, most notably the 12 
spins (six equivalent nuclei and six equivalent ^^C nuclei) in ^^C labelled 
benzene |132fl33fl34] . A partial review of experiments in this field has now 
been published |135] , and there has also been significant theoretical work |136j . 
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2.6 Quadrupolar nuclei 



Nuclear spins with spin quantum numbers greater than one half have quadrupo- 
lar interactions with electric field gradients |137|138j in addition to the Zee- 
man interaction with magnetic fields and dipolar and scalar couplings, and 
so are called quadrupolar nuclei. In the liquid state these terms are averaged 
by molecular motion, although their effects are frequently seen in the short 
relaxation times of many quadrupolar nuclei. In solid state samples quadrupo- 
lar interactions can be comparable to or even larger than Zeeman interactions 
and have huge effects on the NMR spectra. As before, the use of liquid crystal 
solvents provides a convenient middle way, allowing the quadrupolar interac- 
tion to be introduced in a measured and controllable fashion, although short 
relaxation times remain a problem. 



More fundamentally still, such nuclei have more than two spin states and so 
do not provide a natural implementation of a qubit (although the behaviour of 
direct NMR observables can often be described using a Bloch sphere picture, 
more detailed consideration shows that the underlying behaviour is consid- 
erably more complicated). Two approaches are then possible: either the spin 
can be used to represent a larger quantum information carrier (for example, a 
spin-1 nucleus such as is a quantum three-level system and so can be used 
to implement a qutrit |139] ). or a group of levels can be used to represent two 
or more quantum bits. This second approach has dominated, with the four lev- 
els of the spin- 1 nuclei ^SNa 11400141014201430144014501460147] and ^Li [TiS] 
being used to represent two qubits, and the eight levels in the spin-| nucleus 
^^■^Cs being used to represent three qubits |149fl5O0151fl52] . These experi- 
ments have used liquid crystal solvents, but studies have also been performed 
on ^^Na nuclei in a single crystal of NaNOa |153] . 



A major disadvantage of approaches of this kind is that embedding multiple 
qubits in a single spin can make it very hard to scale up the system. In 
conventional quantum computing another qubit can be added to the system 
simply by adding another spin-| particle which is coupled to the existing spin 
system. By contrast, to move from two qubits to three requires replacing a 
spin- 1 nucleus with a spin-| nucleus, while a four-qubit system would require a 
(non-existent) spin- ^ nucleus. This weakness can to some extent be overcome 
by using molecular magnets, which use clusters of spins to create systems 
equivalent to very high spin nuclei |154fl55] . but this does not really address 
the fundamental problem. 
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2. 7 Electron spms: ESR and ENDOR 



There have been many proposals for implementing quantum information pro- 
cessing with electrons, or electrons and nuclei, in solid state systems p-56j . 
Here I will only consider a small number of systems which have been studied 
with relatively conventional ESR techniques |157j . 

One system which has been extensively studied |158|159fl60|161|162|163|164] 
is the "^83/2 electronic state of nitrogen endrohedrally doped into a Ceo fullerene 
cage (i-NCeo, or NQCeo), which provides a spin-| system coupled to either a 
spin-| ^^N nucleus or a spin-1 ^'^N nucleus. However, although several quantum 
phenomena have been demonstrated this system has not yet been used for 
implementing full quantum computation. 

An alternative spin system, which should perhaps be better suited to simple 
quantum computations, is provided by the single electron {S = |) and ^H 
nuclear spin (/ = i) in the 'CH radical formed by x-ray irradiation of a single 
crystal of malonic acid |165fl66] . So far, however, this system has only been 
used to demonstrate two qubit gates used to prepare Bell states |167] . ENDOR 
techniques have also been used to demonstrate the spin-bus or 5'-bus proposal 
for quantum computing in CaF2 doped with Ce^"*" [168] and the transfer of 
quantum states between electron spin and ^^P nuclear spin states for use as a 
quantum memory |169j . 



2.8 Other quantum technologies 



There is a very wide range of other technologies which have been suggested 
for implementing quantum computing, but it is not sensible to attempt to 
summarise these here. An excellent recent set of reviews has been published 
describing several major approaches, including trapped ions |17Uj . ultracold 
atoms trapped in optical lattices |171] . superconducting qubits |172j . spins in 
exotic systems such as semiconductors |173j . and possible approaches for a 
quantum internet jl74] . 



3 Logic gates 



Just as classical computers are built out of a small set of basic logic elements, 
quantum computers can be constructed from a universal set of quantum logic 
gates |ll75j . A traditional set of gates for discussing classical computation is 
AND, OR and NOT, as any desired binary logic function can be assembled from 
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these |176] . In fact this set is excessive as all three of these gates can themselves 
be constructed as circuits of nand or NOR gates, and so both nand and NOR 
are universal for classical computation. These two gates are also relatively 
simple to implement in electronic circuits. 



These gates are, however, completely unsuitable for quantum computations, 
which proceed by a series of unitary transformations of quantum states. Uni- 
tary transformations must be logically reversible (since they have inverses), 
and both nand and NOR are logically irreversible, as it is not possible to re- 
construct the inputs knowing only the outputs. It is possible to implement 
classical computations entirely reversibly |176] . and the three-bit toffoli 
gate, or controlled-controUed-NOT gate, is a reversible equivalent to nand, 
and is universal for classical reversible computation. It is not, however, uni- 
versal for quantum computation, as it does not interconvert eigenstates and 
superposition states of qubits, and so confines qubits to a classical subspace. 

Access to superposition states appears initially to create a problem for describ- 
ing a universal set of quantum logic gates, as there are an infinite number of 
superposition states of the general form Eq. [H and it is not possible to create 
an infinite number of different superposition states using only finite circuits 
built from a finite number of logic gates. It is, however possible to find a circuit 
which approaches arbitrarily close to any desired unitary transformation using 
a number of logic gates which rises only modestly with the desired accuracy 
|177|178|179] . Indeed it can be shown that almost any two-qubit gate is uni- 
versal in this pragmatic sense }180|181|182j . In practice, however, it is usually 
more convenient to enlarge the set of gates, and one common universal set is 
made up of three gates. Two single-qubit gates, the Hadamard gate (H) 



— H — 




(19) 



which performs the transformations 

|0)^ 



-) ^ |0) 



and the fourth root of Z gate (T) 



|1) 



(20) 
(21) 



— T — 




(22) 



can be combined to make any desired single-qubit gate. With the addition of 
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the two-qubit controlled-NOT gate (cnot) 




10 



(23) 
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any desired gate can then be constructed with relative ease |175II178] . 



Important single-qubit gates include the three elementary gates X, Y, and Z, 
corresponding to the three Pauli matrices cr^, (Jy, and a^. Of these three gates 
only X is a classical logic gate, being simply a not gate, while Y and Z only 
make sense on a quantum computer. As discussed below, however, it is usually 
more sensible in NMR to implement X and Y gates (and any other rotations 
in the xy-plane) using RF pulses, and then construct H and T from networks 
of these pulses. 

3.1 Global phases and logic gates 

Just as global phases can be ignored when considering qubit states, they can 
also be ignored when designing quantum logic gates, as the effect of a global 
phase in a logic gate is simply to impose an irrelevant global phase on any state. 
Another way of looking at this is that a quantum logic gate is a propagator 
U, that will be implemented by applying some effective Hamiltonian for a 
time t, such that U = exp(— i'H/f:). A global phase in U then corresponds to a 
constant offset in "H, in effect moving the zero of the energy scale. 

The quantum logic gates used in NMR descriptions usually differ from their 
traditional forms in precisely this way, because most NMR treatments place 
the energy zero between the two spin states, treating the Zeeman effect and 
scalar couplings as splittings around this centre, while most other implementa- 
tions of quantum computing place the energy zero coincident with the lowest 
eigenstate. Thus, for example, the fourth root of Z gate in NMR implementa- 
tions typically takes the form 



rather than the form above, Eq. [22l In single qubit systems this is obviously 
completely irrelevant, and it is also irrelevant when a single qubit gate is 
applied in a system with two (or more) qubits, as the gate then takes the form 




(24) 



18 



T (8) 1, if applied to the first qubit, or 1 (8) T if applied to the second qubit. In 
either case the global phase remains a global phase. 



With controlled gates, however, it is necessary to be much more careful. Con- 
sider for example the controUed-T gate, which traditionally takes the form 



— T — 
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(25) 



It might seem tempting to replace this by an NMR-like form, with a phase 
shift on the T-gate, 
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(26) 
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\0 

but this is entirely wrong, as the phase shift is now a local phase shift, and 
these two gates are quite different. It is, of course, permissable to use an 
alternative controUed-T which differs only by a global phase, such as 
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(27) 



as long as the global phase shift is done correctly. 



3.2 Basic methods: single qubit gates 



As noted above, it is only necessary to implement single-qubit gates, which 
only change the state of a single qubit, and one non-trivial two-qubit gate. 
Here non-trivial means that the final state of at least one of the two qubits 
involved depends on the initial states of both qubits, so that the two-qubit 
gate encodes some sort of conditional logic. Non-trivial two qubit gates can 
interconvert product states and entangled states, and so are essential to access 
the full range of possible states. Trivial two qubit gates can either be written 
as products of single qubit gates, or as products of single qubit gates and the 
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(trivial) two-qubit SWAP gate, which simply swaps the states of two qubits 
PIT83] . 



Single-qubit gates correspond to rotating a single spin in its own one-spin 
Hilbert space. For a one qubit computer, implemented using a single nuclear 
spin, this can be achieved by applying RF fields. A completely general rotation 
can be described by three angles 

U {9, ip, (p) = exp[— i6'(/^. sin ip cos (p + ly sin ipsincp + 1^ cos ip)] (28) 

where ip describes the co-latitude of the rotation axis and product operator 
notation |116] has been used. For simplicity it is often best to consider only 
resonant RF fields, so that = it/2 and the rotation is reduced to the con- 
ventional form 

U{e,(f)) = exp[-ie{I^ cos (f) + ly sin (f))] (29) 
where 6 and (p are the pulse nutation and phase angles. 

Pulses of this kind provide a direct route to the basic gates X and Y, which 
can be implemented as 180° and 180° rotations respectively. Rotations about 
axes not in the xy plane can then be implemented as sequences of pulses. 
Rotations around the z axis are easily constructed with composite z-rotations 
|184] . using the identity 

9, = 90; 90° ^ (30) 

where the pulse sequence is written with time running from left to right, so 
that the leftmost pulse is the first pulse applied. This is, of course, the reverse 
of the order used to describe a sequence of propagators 

exp(— 16'/^) = exp(i7r/2 Jj^) x exp(— i6'/a;) x exp(— i 7r/2 J^) (31) 

which can lead to considerable confusion, as it is not always obvious whether a 
particular equation is describing a sequence of propagators or of pulses. Some 
authors have tried to sidestep this by describing inverse propagators, as these 
multiply in the same order as pulses are listed, but this can simply deepen the 
confusion. 

In larger spin systems it is necessary to implement logic gates in a qubit- 
selective manner, so that only a single-qubit is affected by the RF field. In a 
fully heteronuclear spin system qubit selection is simple, as every spin will be a 
long way from resonance with every other spin, and simple hard pulses applied 
on resonance can be used [T3], but in homonuclear systems it is essential to 
be more careful. 

First, it is necessary to ensure that only one spin is directly affected by the 
pulse, which can be achieved by the use of shaped selective pulses p^ll3fl85] , 
with Gaussian |186j , Hermite |187j and BURP |188] pulses being particularly 
popular. Although many more sophisticated selective pulses are known, these 
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must be used with great they are normally optimised for particular 

transformations on the Bloch sphere (such as excitation or refocussing), while 
quantum logic gates must be implemented using general rotations, which act 
properly on all possible initial states. Alternatively, it is possible to use vari- 
ations on "jump and return" and related pulse sequences |189|190j to achieve 




selection |M|74|191|192j . but this is normally possible only when there are just 
two spins of a particular species. 

Secondly, it is also important to think about the free evolution of any unexcited 
nuclei in the spin system. Even if an apparently ideal selective pulse is used, 
any other spins will still evolve during the pulse duration at their own Larmor 
frequencies, picking up phases reflecting both Zeeman interactions and spin- 
spin couplings. With real pulses these spins will also experience additional 
rotations from transient Bloch-Siegert shifts |193j . and all these phases must 
be considered. Early papers only considered the effects of the direct Zeeman 
interaction, neglecting the two smaller terms, but correction of these (not 
always small) terms has become a major topic in more recent work [29|3T] . 




The simplest approach to correcting the direct Zeeman interaction is to choose 
the length of the selective pulses so that the phase acquired is a multiple of 
271, allowing it to be ignored [12]. Again, this approach can be used only when 
there are just two spins of a particular nuclear type as the pulse sequence is 
synchronised with the stroboscopic alignments between the two spins occur- 
ring at time intervals of 1/ Su, where is the difference in Larmor frequencies. 
Methods based on jump and return sequences implicitly use this stroboscopic 
effect, and so can include the desired corrections automatically. 

With more than two spins a more sophisticated approach is necessary, and 
the most common is the use of abstract reference frames [lllj . Rather than 
using a single rotating frame for each nuclear species, a separate frame is 
used for each individual spin, with every spin being in resonance with its 
own frame. It is not, of course, necessary to have a separate RF transmitter 
for each spin, as "virtual transmitters" can be created using phase ramps 
[185] to shift the transmitter frequency and calculating the appropriate initial 
phase for each spin at each point in time. Note that when applying the same 
single qubit gate to two or more spins it is no longer possible to use the 
simple direct method of applying a hard pulse; instead it is necessary to use 
a sequence of selective pulses, exciting the spins one at a time, or a specially 
designed selective pulse addressing multiple spins. Abstract reference frames 
can also be used to correct transient Bloch-Siegert shifts, by rotating each 
frame appropriately to absorb the additional phase shifts. This process is 
conceptually straightforward but computationally tedious, and so is frequently 
handled by a pulse sequence compiler as described below. Undesired evolution 
under spin-spin couplings can be handled in a similar way, although this 
process is slightly more complicated |194] . 
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Abstract reference frames simplify the implementation of single qubit gates by 
avoiding the necessity of correcting phase shifts (^-rotations), but also provide 
a simple way of implementing desired 2;-rotations. Rotating the frame instead 
of the spin is simply a matter of computational book keeping, and so can in 
principle be done instantaneously and perfectly. For this reason it is often sen- 
sible to decompose single-qubit gates into circuits involving 2;-rotations wher- 
ever possible. One important example is the Hadamard gate, Eq. [191 which can 
be implemented using the pulse sequence 90_y ISO^ or equivalently ISO^ 90^^. 
The z-rotations can be absorbed into the reference frame, and so Hadamard 
gates can be replaced by h = 90^, sometimes called the pseudo-Hadamard 
gate [T2p9] . or its inverse, h^^ = 90-y. As Hadamard gates frequently occur 
in pairs, these pairs can be replaced by one h~^ and one h gate, with the ISO^ 
rotations cancelling out; when they do not occur in pairs it is still possible to 
replace them with pseudo-Hadamard gates as long as the ISO^ rotations are 
absorbed into the reference frames. 

An alternative approach which has grown in popularity is to use methods from 
optimal control theory to develop shaped pulses which both provide selective 
excitation and refocus undesirable interactions. This approach is explored in 
more detail in Section 13.71 below. 



3.3 Basic methods: controlled gates 



Next I turn to non-trivial two-qubit gates, and will begin by considering two- 
spin (two-qubit) systems. It is in principle only necessary to construct a single 
gate of this kind, and the traditional gate used in most theoretical descriptions 
is the controlled-NOT gate, or controUed-X gate, Eq. [221 In NMR experiments, 
however, the key two-qubit gate is the controUed-Z gate 
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(32) 



which is easily converted to a controUed-NOT gate 



— X — 



— H 



(33) 
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with the application of a pair of Hadamard gates [21]. Note that, just like pulse 
sequences, gate networks are written with time running from left to right. 



22 



Unlike the controlled-NOT gate, controlled- Z is symmetric between the two 
spins, which is the form expected for a spin-spin interaction, and it can be 
easily decomposed with product operators |116] 



controlled-Z = exp[— i (7r/2) ( 
= exp[-i(7r/2)( 



-lE + h + S,-2hS,)] 



(34) 
(35) 



where the choice between the two decompositions is simply a matter of conve- 
nience. All four terms commute, and so can be considered individually. The 
term is just a global phase, and can be ignored as usual. Terms in and Sz are 
just single qubit rotations, and can be implemented with single-qubit gates, 
or simply absorbed into the reference frame. This leaves the only important 
term, which in Eq. [34] corresponds to evolution under the (Ising like) spin-spin 
coupling term, 7iJ2IzSz for a time 1/2 J; note that there is a close analogy 
between controlled logic gates and spin-state-selective coherence transfer se- 



quences |195fl96] in conventional NMR. The spin Hamiltonian will include 




both Zeeman and coupling terms, but conventional spin-echo sequences can 
be used to remove the undesirable terms [70]. Similar methods can, of course, 
be used to generate closely related gates [T^PT] . such as the controlled-T gate. 



Similar methods can also be used in systems with more than two spins; the 
key element lies in the use of spin-echo sequences to isolate the desired spin- 
spin coupling, while suppressing all undesired interactions. Conventional ap- 
proaches, based on nesting spin-echoes, become very complicated in large 
fully-coupled spin systems, and this should instead be done efficiently using 
sequences based on Hadamard matrices |110fl97fl98] . In a system of n spins, 
all of which are coupled weakly to one another, any desired coupling can be 
isolated with around pulses, although it is not always necessary to refocus 
all the additional couplings |194] . 

Most spin systems are not fully coupled, and a much smaller number of pulses 
can be used to isolate desired couplings in such systems |110fl97j . If a desired 
coupling does not exist, or is too small to be useful, then it must be generated 
indirectly. This is straightforward in principle, as quantum swap gates [15|183] 
can be used to move logical qubits around a spin system to bring them into 
contact; in practice it is not necessary to use full SWAP gates, and considerably 
more efficient approaches are known |82|199j . It is, of course, essential that 
every pair of spins must be connected either directly or indirectly by some 
chain of usable couplings: it is not possible to perform two-qubit gates between 
spins in completely separate spin systems. 

More complex gates, such as the toffoli gate, or controUed-controlled-NOT, 
can then be constructed using quantum circuits |175fl78j . but it is often better 
to construct these gates directly from NMR primitives, as this can lead to 
significantly simpler implementations |2UUf2Ul] . 
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3-4 Transition selective pulses 



An alternative, conceptually simpler, approach to implementing controlled 
gates is simply to use transition selective pulses |202j . The NMR spectrum of 
a two-spin system will contain two doublets, and a sufficiently long selective 
pulse will only excite one transition in a doublet, and so will only affect the 
state of the excited spin if the coupled partner is in the desired state. It 
might seem that applying a 180° pulse to one transition would implement a 
controlled-NOT gate, but it is necessary to think carefully about global phases. 
The propagator corresponding to a 180° rotation is —icTx, which differs from 
a NOT gate by a global phase, reflecting spinor behaviour [9] . For a transition 
selective pulse \203\ this phase is no longer global, and so cannot be ignored: 
the propagator is 
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which differs from a controlled-NOT gate, Eq. [221 by a ^-rotation on the control 
spin. 

In a two-spin system this approach has no particular advantage, and is indeed 
fundamentally equivalent to the more traditional approach based on pulses and 
delays, which can be considered as a "jump and return" style implementation 
|192] of a transition selective pulse. Alternatively, the transition selective pulse 
method relates to the pulses and delays method much as the selective popula- 
tion transfer experiment |204j relates to INEPT [205j. In multi-spin systems, 
however, transition selective pulses provide a simple and direct way of im- 
plementing more complex quantum logic gates. For example, in a three-spin 
system the toffoli gate can be implemented directly as a transition selective 
pulse addressing one of the four lines in the double doublet corresponding to 
the target spin |^48|69|71|72j . A corresponding disadvantage of this method 
is, of course, that constructing a two qubit gate in this system will require a 
pulse which selects two of the four transitions for the target spin [21] . 

3.5 Composite pulses 

Composite pulses |206j have found widespread use in conventional NMR ex- 
periments to reduce the effects of a wide range of experimental imperfections, 
most notably off-resonance effects and pulse length errors arising from RF 
inhomogeneity. Similar imperfections are likely to affect most experimental 
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implementations of quantum information processing and there has been con- 
siderable interest in applying these ideas. Most conventional composite pulses 
are not, however, suitable for use in quantum computers, as they are opti- 
mised for particular classes of initial state: for example, most composite 180° 
pulses are optimised either for inverting the population of a spin, or for pro- 
ducing a spin echo on a coherent superposition. By contrast, pulses used on 
quantum computers must be general rotors, which perform well for any initial 
state. Composite pulses of this kind are rarely used in conventional NMR, 
but a small number of so-called Class A composite pulses [206j are known, 
and these have been developed for wider use PT|207II2U8] . A method for con- 
structing general rotors from conventional point-to-point pulses has also been 
described [209] . 

An early composite 90° pulse tackling off-resonance errors was described by 
Tycko |210j . replacing a 90° pulse with the three pulse sequence 385°320!_2,25°, 



has been generalised to give the CORPSE family of composite pulses |208II211] , 

in which a 6x pulse is replaced by three pulses, applied along the +x, —x and 
+x axes as before, with flip angles given by 



where rii, n2 and 71,3 are integers, with the best results [208] occurring for 
rii = n2 = 1 and 77.3 = 0. These sequences have been demonstrated by NMR 
[211] and SQUID [212] experiments. Pulse sequences have also been designed 
which are tailored for particular off- resonance effects [191] . 

There has been considerably more interest in composite pulses to tackle pulse 
length errors, which can largely be traced back to a three pulse composite 
180° pulse due to Tycko and coworkers [210] or to the BBl family of sequences 
discovered by Wimperis [213] . Tycko's pulse sequence has been generalised to 
give the SCROFULOUS family of composite pulses [208] . but using the BBl 
family is preferable in most cases. 

BBl differs from many other composite pulses in that it seeks to design an 
error-correcting pulse, which can be combined with the naive error-prone pulse 
to give a more accurate compound pulse, much as a contact lens can be used 
to correct eyesight. Originally [213J this error correcting sequence (sometimes 
called a Wl sequence) was placed before the naive pulse, but it can instead 
be placed after the naive pulse, or indeed in the middle of it [208f214] . It 
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comprises three pulses, in the form 180^^ 360^^ ISO^i) with 02 = 3 0i and 

01 = ± arccos 

where the choice of sign is unimportant as long as it is made consistently. 
The method can be extended to build two-qubit gates which are robust to 
variations in the size of the underlying scalar coupling [215,216,214] 

BBl has proved a remarkably successful composite pulse, and is surprisingly 
difficult to improve upon. The quality of a composite pulse for quantum com- 
puting can be assessed in various ways, but in practice there are two important 
families of approaches. The most direct approach is to expand the propagator 
for the composite pulse as a power series in the size of the error, and deter- 
mine the size and order of the lowest order error term. As an example consider 
implementing a 180° pulse using a naive pulse with a fractional pulse length 
error of e, so that the flip angle of the pulse is in fact 180° x (1 + e). The 
propagator is then 

/ -i\ /7r/2 \ 
exp[-i7r(l + e)/J= + ) (^1) 

1^-i ; [O 7r/2) 

and so the naive pulse has an error of order e. Alternatively the quality can 
be assessed by calculating the propagator fidelity F between the desired prop- 
agator U and the actual propagator V, given by 

F =\Tt{VU-'^)/Tt{UU~'^)\, (42) 

and then expanding the fidelity as a power series in the error. For the naive 
pulse considered above the fidelity is 

F=l- e^TT^S + 0{e^) (43) 

and so the naive pulse has infidelity of order e^. The difference between these 
two methods of assessing a pulse must be borne in mind when comparing 
pulses in different papers; in general an error of order n will correspond to an 
infidelity of order 2n. 

BBl pulses can be derived by designing composite pulses which suppress first 
order pulse length errors, but it turns out that BBl also suppresses second or- 
der errors automatically, leaving only third order errors (sixth order infidelity). 
It is not clear why this fortuitous double cancellation occurs; it is not a general 
feature of composite pulses. Other pulses with similar properties are known 
|217j . but these have no advantages over BBl. Beyond this, BBl pulses are 
also relatively robust to off resonance-errors |208j , and generally insensitive to 
small errors in their implementation, so that BBl pulses work in practice very 
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much as expected from theory [214j . In addition to NMR experiments [6T|214] 
BBl pulses have been demonstrated in electron spin resonance |160fl63j and 
have inspired applications in other fields |218f219f22Uf221] . 

Although BBl has proved highly successful, it is obviously interesting to seek 
still better pulse sequences, and Brown et al. have tackled this in two ways 
[222f223] . Firstly they have shown how the BBl approach can, in effect, be 
nested, creating ever higher orders of simultaneous correction. A robust 90° 
pulse from the B4 family of pulses (which remove the third order error term) 
has been implemented in NMR experiments |214j , but this composite pulse is 
very long (the correction sequence contains 27 pulses with a total length equiv- 
alent to a 7200° rotation) and does not perform much better than BBl. Sec- 
ondly they have described a general method, using insights from the Solovay- 
Kitaev theorem |179j . showing how arbitrarily accurate composite pulses can 
be constructed in general by building a series of correction sequences which 
correct errors one order at a time. An expanded version of part of their method 
written in more conventional NMR notation is also available |224j . These ideas 
have subsequently been extended to multi-qubit systems [225] . 



3.6 Strongly modulated composite pulses 

Strongly modulated composite pulses [89|102f226j are not composite pulses in 
the conventional sense, but are in many ways closer to shaped pulses [185] . 
Conventional composite pulses are comprised of a small number of pulses 
which usually vary only in length and phase, while shaped pulses are made up 
from a large number of pulses of the same length, with control of both ampli- 
tude and phase. Strongly modulated composite pulses contain a small number 
of pulses, but these pulses can differ in amplitude and frequency, as well as 
length and phase. In practice strongly modulated pulses are always imple- 
mented as shaped pulses, with frequency shifts implemented as phase ramps 
[185] . so that strong modulation can be thought of as an unusual method for 
parameterising a shaped pulse, and strongly modulated pulses are designed by 
numerical simulation rather than analytical calculation. Unlike most conven- 
tional methods, however, the propagator for the pulse can be calculated using 
the short composite pulse description, rather than the much longer shaped 
pulse description, simplifying calculations. It is, of course, necessary to ensure 
that the sufficient digitisation is used when implementing the composite pulse 
as a shaped pulse. 

Strongly modulated pulses were originally developed to allow desired quantum 
logic gates to be implemented in complex multi-spin systems. The conventional 
approach, using spin-echoes to isolate desired terms in the spin Hamiltonian 
and then constructing a logic circuit from these elementary terms, can become 
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impractical in larger spin systems, due both to the large number of interac- 
tions which have to be considered and interactions between non-idealities aris- 
ing from individual pulses. Strongly modulated pulses seek to implement the 
desired logic gate in one go, simply optimising the propagator of the whole 
composite pulse, calculated from the sum of the background and RF Hamil- 
tonians, by changing the parameters describing the RF component |89]. This 
also allows the pulses to correct for any non-idealities in the system: for exam- 
ple, strongly modulated pulses can deal with the slight deviations from weak 
coupling found in homonuclear systems, by simply using the full coupling term 
in the Hamiltonian used to calculate the pulse propagator. In the same way 
non-idealities, such as transient Bloch-Siegert shifts and interactions between 
selective pulses, are automatically incorporated in the calculation and cor- 
rected for in the pulse sequence, as can effects such as short delays around 
pulses imposed by the spectrometer hardware. The effects of incoherent ef- 
fects, such as RF inhomogeneity, are not automatically incorporated, but can 
be included by simultaneously optimising the propagator for a range of RF 
strengths [89)1227] : this also allows incoherent effects to be distinguished from 




decoherent effects |226j . 

The discussion above assumes that a strongly modulated pulse can be found 
to implement a desired propagator. As the NMR Hamiltonian supplemented 
by RF fields provides a universal set of quantum controls, it should always be 
possible (neglecting RF inhomogeneity) to find a suitable pulse, as long as a 
sufficient number of pulses and total pulse length are available. The required 
number of pulses is usually found by trial and error, with the number of 
pulses being increased until a suitable solution is found. The total pulse length 
required can be estimated from optimal control theory (see below), and then 
optimised by the algorithm, or can simply be found by trial and error. It is 
usual to start optimisation from a fairly smooth low power initial pulse, and 
to use several initial pulse sequences, retaining those which seem to converge 
to an acceptable solution for further optimisation. If desired penalty functions 
can be added to dissuade the algorithm from designing pulses with excessive 
length, offset frequencies, or RF power, but this is not always necessary. In 
practice it is usually possible to find suitable pulse sequences fairly easily, 
which perform desired transformations with high fidelity (values above 0.999 
are usually possible). When the effects of RF inhomogeneity are included 
the process becomes a little more complicated, but it is often possible to find 
sequences with fidelities above 0.99 over a reasonable range of inhomogeneities. 

Strongly modulated pulses have been demonstrated in a range of NMR quan- 
tum information processing experiments with spin-i nuclei in the hquid state 
[I0]l89ll9ll971l98]ll02y226y227y228y229y230] or solid state p^0llT2T] . as well as to 




strongly coupled systems [128], quadrupolar nuclei [147fl53] and ENDOR 
[231] . Until the arrival of GRAPE techniques, described below, it seemed 
likely that strong modulation would become the dominant method for de- 
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signing pulses for NMR quantum information processing. There is, however, 
some concern as to whether using such pulses is in some sense cheating, as their 
design requires a full simulation of the entire spin system, and the difficulty 
of the simulation increases exponentially with the number of spins involved 
[iO] . Indeed, the best way to design strongly modulated pulses would be to 
use a quantum computer! For this reason some authors have avoided strongly 
modulated pulses where possible, and have selected spin systems carefully to 
minimise these problems [232]. 



3. 7 GRAPE and other methods from optimal control 

As strong modulation is simply an unusual method of parameterising a shaped 
pulse, allowing the numerical optimisation of arbitrary transformations with 
arbitrary Hamiltonians, one might wonder why shaped pulses shouldn't be 
designed by brute force numerical optimisation of the amplitude and phase 
of each period, with no attempt at imposing a particular parameterisation. 
Historically this has not been done, due to a perception that optimising over 
a large number (hundreds or even thousands) of amplitudes and phases would 
be computationally infeasible. 

Despite these perceived difficulties there has always been some interest in this 
approach, arising from the field of optimal control theory. Once it has been 
established whether a transformation is possible in principle |233j . the next 
obvious question is how the required resources, such as the time required to 
implement the transformation, can be minimised. While there has been some 
abstract work on this question from the viewpoint of quantum computation 
|234|235|236j a more profitable approach for pulse design is provided by the 
field of coherent control |237f238f239] . Early applications to NMR largely stud- 
ied methods for performing coherence transfers between spins in the minimum 
possible time |240f241|242] . including relayed transfers via intermediate spins 
|243] : these ideas have obvious applications in branches of NMR spectroscopy 
where signal loss due to spin-spin relaxation during coherence transfers is a 
problem |244f245f246j . 




These ideas were then generalised to the implementation of quantum logic 
gates |247j and the grape (gradient ascent pulse engineering) algorithm was 
described [123] which allows arbitrary shaped pulses to be developed to per- 
form desired unitary transformations by adding control fields to a background 
Hamiltonian. The key idea behind grape is a simple and efficient method 
[123] for estimating the gradient of the fidelity between a shaped pulse and 
the desired transformation, thus allowing conventional optimisation methods 
to be used to design pulses. The traditional method for estimating a gradient 
for a composite pulse with n independent time periods, and so described by 
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2n variables, requires 2n + 1 complete evaluations of the pulse propagator, 
but the GRAPE approach enables the gradient to be estimated using only 2 
complete evaluations, by eliminating unnecessarily repeated calculations. This 
makes it possible to optimise shaped pulses described by very large numbers 
of parameters (hundreds, or even thousands of time periods). 

GRAPE methods have been used to design a wide variety of exotic shaped 
pulses with applications in conventional NMR |248')249I2 50|251] and to quan- 
tum information processing using both NMR p^22 252 253] and other tech- 
niques |221|231|254] . An excellent recent description of the approach can be 
found in [21]. Two particularly interesting developments are grape pulses 
which work entirely by phase modulation [250], and give compensation for RF 
inhomogeneity over extremely wide ranges, and pulses which work entirely by 
amplitude modulation (except for 180° phase shifts) |231] and so can be used 
in experimental fields where small angle phase shifts are not available. 

Given that grape pulses are in principle more powerful than strongly mod- 
ulated composite pulses, one might expect them to become the method of 
choice in NMR quantum information processing, and there are some signs of 
this happening. However, grape pulses are vulnerable to the same fundamen- 
tal criticism as strongly modulated composite pulses, namely that their design 
requires a full simulation of the entire spin system. An attempt has been made 
to mitigate this problem by combining grape methods with pulse sequence 
compilers, as described below. 

3.8 Dynamical decoupling 

Dynamical decoupling |255|256|257j extends the ideas of coherent control to 
the problem of controlling open quantum systems, such as quantum systems 
undergoing decoherence. The ideas can be related to ideas in conventional 
NMR |257f258j . including the use of decoupling techniques to suppress un- 
wanted couplings, and the use of spin locking to create an effective relaxation 
time Tip. These ideas have been extensively studied |259f260f261f262f263j and 
have found applications in NMR [264j and EPR |162] . Conceptual links can 
also be drawn between dynamical decoupling and the use of decoherence free 
subspaces, discussed in Section [5751 below. 

3.9 Pulse sequence compilers 

The pulse sequences used in NMR quantum information processing are, even 
by the standards of NMR, extremely long and complicated. For example, one 
pulse sequence used in an NMR implementation of Shor's quantum factor- 
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ing algorithm contains 299 shaped pulses [113] . While some solid state NMR 
experiments can involve larger numbers of pulses, these are usually found 
in highly repetitive decoupling or mixing sequences, rather than individually 
placed pulses. The sheer complexity of writing such pulse sequences has led 
many groups to develop simple computerised tools, usually called pulse se- 
quence compilers, which can handle many of the more mundane steps in the 
calculation. 

Pulse sequence compilers can act at a variety of levels, from the abstract net- 
work description, where they can be used to cancel unnecessary gates |113j . to 
low level frame tracking, where they are used to keep a record of frame rota- 
tions and extraneous couplings PT|194] . The combined use of GRAPE pulses, 
designed to work within a spin subspace, and pulse compilers which keep track 
of interactions with spins outside the subspace seems particularly powerful, 
and is described in detail elsewhere [3T| . 



3.10 Geometric phase gates 

Geometric phases are a topic of considerable interest in many branches of 
physics |265 '.'^266,267,268,269]. and have over many years been explored in the 
context of NMR ^270.271.27211775] , ESR [2711775] and NQR |276M277a278M279M280M28H282Eg5] . 
The simplest example of a geometric phase is Berry's phase |266j . which arises 
when a Hamiltonian is adiabatically varied around a circuit. A quantum sys- 
tem which starts in an eigenstate of the initial Hamiltonian will always remain 
in the corresponding instantaneous eigenstate of the time- varying Hamiltonian 
(as the process is adiabatic) and so will return to its initial state (as the pro- 
cess is cyclic). The state can, of course, acquire a phase in this process, and 
as well as the conventional dynamic phase, which depends on the average en- 
ergy of the state and the time taken for the process, the state will acquire an 
additional geometric phase, whose size depends only on the solid angle sub- 
tended by the circuit in some appropriate parameter space, and in particular 
is independent of the time taken to complete the circuit. 

A simple way of thinking about these geometric phases is that the motion 
of the Hamiltonian leads to a fictitious magnetic field (a gauge field) |276] 
which interacts with the spins. The size of the field is proportional to the rate 
of change of the Hamiltonian, and so the product of the field strength and 
the cycle time is constant. This picture makes clear that different spin states 
acquire different geometric phases, so that the Berry phase is not a global 
phase, and the difference between two geometric phases can be detected as a 
relative phase shift in the evolution of a superposition. Alternatively, geometric 
phases can be detected in interference experiments; in effect these work by 
replacing the global phase by a controlled phase, controlled by the "which 
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way" path information |271] . 



Berry's phase can be generahsed in several important ways. In NQR experi- 
ments the degeneracy of spin states means that the system need not return 
to the same state at the end of a cychc evolution; instead mixing can occur 
within the degenerate spin states. Unlike phase shifts, this more general mix- 
ing need not be Abelian |277] . which can lead to more complex behaviour |278] 
than that seen for spin-| nuclei. More usefully, Aharonov and Anandan have 
shown [284] that Berry's phase can be seen as the adiabatic limit of a much 
more general phase, and that in non-adiabatic cases this phase depends on the 
trajectory of the state, not the trajectory of the Hamiltonian. Such Aharonov- 
Anandan phases are relatively simple to observe [271j : indeed they frequently 
arise in quite conventional NMR experiments. For example, the fact that evo- 
lution under two 180° pulses with phases differing by is equivalent (up to a 
global phase) to a ^-rotation through an angle of 20 can be understood as an 
Aharonov-Anandan phase [27lj. 

In addition to their fundamental interest, it has been widely suggested that 
geometric phases could be technically important, as their size depends only 
on the geometry of the cyclic path and not on other details, and thus geo- 
metric phases might be robust to variations in experimental parameters. It 
is, however, not clear that such robustness can always be achieved. Firstly, 
the geometric phase usually (but not always [285] ) occurs on top of a (usually 
much larger) dynamic phase, and it is necessary to remove this dynamic phase 
term. This can usually be achieved by refocussing it with a spin-echo sequence, 
but the quality of the geometric phase will then depend on the quality of the 
refocussing step. Secondly, the geometry of the path can itself depend on ex- 
perimental parameters, and it is necessary to think carefully about exactly 
what determines the geometry in a particular physical situation. In particu- 
lar, with Aharonov-Anandan phases it is important to ensure that the state 
evolution does in fact form a closed cycle; if it does not then the path must be 
closed with a geodesic, changing the solid angle subtended by the path from 
its naive value. 



Despite these caveats there has been considerable interest in the use of geomet- 
ric phases to implement quantum logic gates. The first experimental demon- 
stration [58|286j used the controUed-acquisition of Berry phases to implement 
a two-qubit quantum logic gate in the heteronuclear two-spin system pro- 
vided by ^'^C labelled chloroform. This has been followed by many experiments 
and proposed experiments in NMR |152|287|288II289|29U|291II292|293] and ESR 

later papers are largely based on an earlier proposal for holonomic quantum 
computation [301] . which uses non- Abelian phase mixing effects). There have 
also been theoretical studies of the robustness of such gates against decoher- 
ence |302] . variations in control parameters |303ll304^305),3TIB] . and non-cyclic 
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evolution [307]. 



3.11 Quantum cellular automata 

The discussion so far has assumed that it is necessary to selectively address 
single qubits and pairs of qubits in order to implement general quantum com- 
putations. In fact it has been known for many years that quantum information 
processing can be implemented in systems with much less control of individ- 

to as quantum cellular automata, or QCA, although that this term is also 
used to describe classical cellular automata implemented with quantum dots 
or other quantum technologies |319|320|321|322|323|324|325] . and, to add fur- 
ther confusion, there is occasional discussion of implementing quantum in- 
formation processing on quantum dot cellular automata, sometimes called 
coherent quantum dot cellular automata |326] . In this discussion I will only 
consider schemes capable of implementing quantum information processing 
without selective qubit addressing; these ideas may have important applica- 
tions to quantum information processing in optical lattices [327]. 

Quantum cellular automata proceed by dividing up spins into two or more 
groups (usually labelled A, B and so on), such that all the spins of a particu- 
lar type can be addressed simultaneously but it is not possible to separately 
address individual spins within a given type. The spins are assumed to interact 
through some network of couplings, typically acting between nearest neigh- 
bours in a linear chain or two dimensional array. The simplest practical scheme 
[M|318] assumes a linear chain of spins, all of the same type and connected 
by Ising couplings. In systems of this kind it is always possible to address the 
end spins separately from the spins inside the chain, and this provides suffi- 
cient control to implement full quantum computing [51]; it is not necessary to 
be able to distinguish the two end spins from each other, although this does 
help. It is even theoretically possible to use a completely symmetric ring of 
spins |317] . but it is not clear that this scheme can in fact be implemented in 
physical systems. 

The simple linear chain scheme has been implemented on an NMR quantum 
cellular automaton using three ^^C nuclear spins in alanine [91], which was 
used to implement both a two-qubit and a three-qubit algorithm. The system 
was also used to explore quantum mirroring, in which qubits are moved along a 
chain to their mirror-image locations. As mirroring moves a qubit from one end 
of a chain to the other it can be used to transfer spin states along a quantum 
wire, thus moving a qubit around a quantum computer |328|329|330|331] . 




In addition to quantum mirroring, there has also been interest in implement- 
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ing particular quantum information processing tasks which do not require 
complete control of a quantum system; indeed, in some cases the implemen- 
tation may be simpler if the spin system has a QCA structure. One obvious 
and important example of this is entanglement assisted spin state measure- 
ment |228|332] . and the related technique of entanglement assisted magnetic 
field sensing |114|115] . where the use of a star-topology QCA permits a highly 
entangled state to be prepared with relative ease. 

3.12 Quadrupolar nuclei and strongly coupled systems 

The methods used to implement quantum logic gates in quadrupolar systems 
are not fundamentally different from those used for spin-| nuclei, but there is 
a difference of emphasis. The large size of quadrupolar couplings means that 
transition-selective pulses are simple to apply, and provide the most natu- 
ral way of implementing controUed-NOT gates |14U|142II148] : they can also be 
used to implement controlled phase-shift gates using geometric phases |152] . 
Much is often made of the relative simplicity of implementing two-qubit gates, 
especially the SWAP gate [148] in such systems, but it must be remembered 
that implementing single-qubit gates can be correspondingly difficult, requir- 
ing pulses which simultaneously or sequentially select multiple transitions. It 
is also necessary to consider carefully the effects of background evolution un- 
der the quadrupolar coupling on spin states not affected by selective pulses 

m- 

The methods used to implement quantum logic gates in strongly coupled sys- 
tems are frequently complicated and depend strongly on the details of the 
system. In spin-systems where the weak coupling approximation is starting 
to break down the use of strongly modulated composite pulses [12011121] or 
GRAPE pulses |122j provides a simple solution. In more strongly coupled sys- 
tems strongly modulated composite pulses |128j have been used, but the use 
of single transition selective pulses |127] or methods adapted from solid state 
NMR |135] may be more appropriate. 



4 Initialisation 

Initialisation is the process of preparing qubits in some well defined initial state 
from which a computation can proceed. This process is not much discussed 
in conventional computation, as it is seen as trivial and uncontroversial, but 
it plays an important and explicit role in reversible computation |176j , as the 
initialisation step is fundamentally irreversible (the bits must be set to the 
desired initial state whatever their state was before the initialisation process). 
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This fact leads to a minimum energy cost for computation [333], resolving the 
apparent paradox of Maxwell's Demon. 

The situation is perhaps even more stark for quantum computers, where Mer- 
min has forcefully argued pi|334] that the state of a qubit is not even a 
meaningful concept until the qubit has been initialised by a measurement. 
Whether or not one accepts Mermin's arguments, it is clear that initialisation 
and measurement are intimately connected. The simplest way of initialising a 
qubit is to measure it in the computational basis: if the result is |0) then the 
qubit should be left alone, while if the result is |1) then it should be flipped 
with a NOT gate. This procedure, however, requires access to a true quantum 
measurement device, which projects the qubit into the measurement basis. As 
we shall see in the discussion on readout, conventional NMR measurements do 
not have this property, and so cannot be used to initialise spin states. Instead 
it is necessary to use indirect methods, in which the spin state is implicitly 
measured by the environment. 



4-1 Preparing pure states by cooling 

The most obvious method for preparing a spin system (or indeed any system) 
into a pure quantum state is by cooling it into its ground state. Considering 
the Boltzmann expression it is immediately clear that this will only work if 
the thermal energy kT is much less than the energy difference AE = hv be- 
tween the ground and first excited state. For NMR with a Larmor frequency 
below 1 GHz this requires T ^ 0.05 K, which is clearly impractical for most 
liquid state samples, although it can be reached by solid state NMR or for 
liquid state studies of ^He [335j . For this reason, most implementations of 
NMR quantum information processing have not used pure states, but rather 
pseudo-pure states as described below. There is also considerable interest in 
developing methods for generating non-Boltzmann populations in spin sys- 
tems [2U], although with the exception of techniques based on para-hydrogen 
[1^50|336|337|l538j these attempts have not yet been completely successful. 
The situation is, of course, more promising for electron spins, due to the much 
higher Larmor frequencies involved, although cooling to the solid state is still 
required [339] . 

One important point that is frequently neglected is that it is not really suf- 
ficient to have just a method for preparing a pure state at the start of a 
computation. Instead, many quantum information processing methods, most 
notably error correction, require the ability to reinitialise at least some qubits 
in the middle of a computation. For this reason, no method yet demonstrated 
is entirely satisfactory, although some methods could in principle be extended 
to permit this. 
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4-2 Pseudo-pure states and entanglement 

Because of the difficulty in preparing pure spin states in NMR systems, almost 
all NMR quantum information processing experiments have used pseudo-pure 
states [71l8|9] . sometimes called effective pure states [10]. There are many dif- 
ferent ways of preparing such states, which will be discussed below, but it is 
useful to begin by considering their general properties. 

A pseudo-pure state in a system of n spins is simply a mixed state of the form 

(l-e)l + e|V^)(^| (44) 

where is the corresponding pure state, 1/2*^ is the maximally mixed state 
(that is, equal probability of the system being found in any particular state), 
and e, the excess probability of finding the system in the desired pure state, can 
be taken as a measure of the purity of the state. Since quantum computations 
proceed through a series of unitary transformations the pure component of a 
pseudo-pure state will evolve in exactly the same way as a pure state would, 
while the maximally mixed state is unaffected by any unitary transformation; 
thus the behaviour of a pseudo-pure state is effectively the same as that of a 
pure state. For initialisation the pure component can be taken to be the 
ground state |0) = |00 . . . 0). 

Significant differences will, of course, be seen at the end of the computation, 
when an attempt is made to measure the final state of the quantum com- 
puter. With conventional quantum measurements the result may correspond 
to some component from the maximally mixed state, rather than the desired 
pure component, and for the low purities found in NMR systems such results 
might be expected to overwhelm the tiny pure component. However, NMR 
measurements are not sensitive to the maximally mixed state (in effect, the 
signals from the different components of the maximally mixed state cancel 
out), and only the signal from the pure component is observed. Of course the 
size of the observable signal does depend on the purity, and so pseudo-pure 
states will give much smaller signals than pure states. 

This small signal size might not immediately concern NMR spectroscopists 
who are accustomed to observing signals from spin systems with very low 
polarisations. It is, however, important to realise that the states observed 
in most conventional NMR spectra are not pseudo-pure states [19], and so 
instinct may not be reliable. This problem was immediately raised by Warren 
[16j, who showed that the maximum pseudo-pure state signal obtainable from 
a system of n identical spins in thermal equilibrium cannot exceed a simple 
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bound IT61I22] of 
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where the approximation apphes in the high temperature limit, that is hu ^ 
kT, which is appropriate in conventional NMR experiments. The key part 
of this resuh is the first term of the high temperature result, which shows 
that the observable signal decreases exponentially with the size of the spin 
system [16]. This effect appears to limit NMR quantum computing based on 
pseudo-pure states to around 10-20 qubits ^lOnlTj . 

This practical concern about the scalability of NMR quantum computers is 
bad enough |,20J, but has also led to more fundamental concerns [IB] as to 
whether NMR quantum computers are really quantum devices at all! When 
considering this claim it is vital to understand that the word quantum is 
being used in a very technical sense, essentially equivalent to provably non- 
classical, to describe systems which transcend known limits on the processing 
power of classical systems, and so the two formulations commonly heard (NMR 
quantum computers are not really quantum, and NMR quantum computers 
are quantum, but not usefully so) are less different than it might at first 
appear. 

The essence of this discussion revolves around the problem of quantifying 
the extent of entanglement in pseudo-pure states. Entanglement is widely 
believed to play a key role in quantum information processing |340j . and it has 
been established that many quantum algorithms must involve the generation 
of entangled states |341j . Unfortunately pseudo-pure states corresponding to 
entangled states (sometimes called pseudo- entangled states) are only genuinely 
entangled if their purity is high enough [IB] . It should be noted, however, that 
although individual NMR spin states can be described using a classical model, 
attempts to describe complete NMR experiments in this way [342] have so far 
proved unsuccessful. 

Problems of this kind were first considered by Werner |343] . and two-qubit 
pseudo-pure states where the pure component is a Bell state are usually known 
as Werner states. These states were analysed in detail by Peres |344j who 
showed that a pseudo-entangled state of this kind is entangled if e > | and 
is separable at or below this bound. (Note that the bound corresponds to 
a fractional population of one-half in the Bell state, not one-third, as the 
two-thirds part of the ensemble in the maximally mixed state can itself be 
decomposed as an equal mixture of the four Bell states.) 

Quantifying the entanglement of pseudo-pure states in larger spin systems is 
not a completely solved problem, but several key results have been derived. 
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Braunstein and coworkers [18] originally derived a lower bound 



12 

e < 46 

below which all pseudo-pure states are certainly separable, and an upper 
bound 

12 
^ ^ 1 + 2"/2 2^ 
above which provably entangled states can be made, leaving an intermediate 
region which is not well understood. Subsequent results have tightened these 
bounds further |345ll346f347j . but they do not as yet coincide. 



4-3 Preparing pseudo-pure states 



Although the basic idea of using pseudo-pure states is straightforward, it is 
necessary to consider how they can be prepared. For a single spin the thermal 
state, with an excess population in the lower level, is a pseudo-pure state, but 
this is not true in larger spin systems [19]. Furthermore, in such systems the 
thermal state cannot be converted into a pseudo-pure state by any sequence of 
pulses and delays, or indeed by any unitary transformation. This is most easily 
seen by noting that the eigenvalues of the two states are quite different, and 
that the eigenvalues of a matrix are invariant under unitary transformations. 
(For a pseudo-pure ground state these eigenvalues are simply the diagonal 
elements of the density matrix, that is the populations of the basis states, 
while a general pseudo-pure state will also be diagonal in some appropriately 
chosen basis.) 

A pseudo-pure state has a single large eigenvalue, e+(l — e)/2"', corresponding 
to the population of the desired ground state, and 2" — 1 identical smaller 
eigenvalues, (1 — e)/2", corresponding to the equal populations of the various 
excited states; by contrast a thermal state has a complex pattern of eigenvalues 
reflecting the varying populations of the different excited states. For example, 
for a homonuclear two-spin system the relative thermal state populations in 
the deviation density matrix are 

{1,0,0,-1} (48) 
while the desired pattern for a pseudo-pure state is 

(49) 

This comparison also makes clear how in principle a pseudo-pure state can be 
prepared: all that is necessary is to average the values of the smaller eigenvalues 
while leaving the largest eigenvalue untouched. 
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4-4 Temporal averaging 



Perhaps the conceptually simplest approach to averaging the smaller eigen- 
values is the method of temporal averaging [55], in which the final result 
is obtained by averaging the experiment over a number of different starting 
states. As quantum logic gates and the NMR detection process are linear in 
the input state, this is entirely equivalent to performing the experiment on a 
single averaged input state. By this means the desired component can be re- 
tained, while the undesirable terms cancel out. Conceptually this is somewhat 
similar to the use of phase cycling, with the averaging occuring over a range 
of different initial states. 

For example, in a two-spin system, a pseudo-pure state can be achieved by 
averaging over the three initial states (described as before by the relative 
populations in their deviation density matrices) 

{1,0,0,-1} {1,0,-1,0} {1,-1,0,0} (50) 

as these average to the desired pseudo-pure state, Eq. HOl This example also 
makes clear which initial states should be averaged: the first state is simply the 
thermal state, while the remaining two can be obtained by cyclically permuting 
the three smaller populations in the thermal state. 

In a system of n spins there will be 2" — 1 smaller populations in the thermal 
state pth, and it clearly suffices to average the 2" — 1 cyclic permutations of 
these populations to obtain the desired pseudo-pure state po- In general one 
may write 



— y: Pjpt^p] (51) 



^ " i=o 

where the Pj are permutation operations, with Pq being the trivial identity 
operation. This approach is called exhaustive averaging, and will work for any 
size of spin system and for any set of thermal populations (so it will work 
for heteronuclear spin systems as well as homonuclear systems), or indeed 
for almost any set of initial populations. All that is necessary is that the 
desired state |0)(0| has the largest population in the system; if this is not the 
case then each permutation must be preceded by an initial operation which 
swaps the largest population to |0)(0|. The desired permutation operations 
can be implemented using gate networks [5511348] . For the two-spin case the 
two permutation networks can be constructed using controUed-NOT gates 



(52) 



and general permutations can always be constructed from not, controlled- 
NOT and TOFFOLI gates, although this theoretical description assumes that 
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permutations can be implemented without errors arising from experimental 
imperfections and relaxation, which may not be realistic in complicated spin 
systems. 

An obvious disadvantage of exhaustive temporal averaging is that the num- 
ber of initial states used increases exponentially with the number of qubits in 
the system. This increases the number of separate experiments by the same 
margin, and thus the overall time taken to implement the calculation; this 
completely cancels any exponential improvements in computational complex- 
ity that can arise from using quantum algorithms. Because of this limitation 
there has been interest in finding more efficient methods, using non-cyclic 
permutations and applying unequal weightings to different permutations. For 
example, preparing a pseudo-pure state in a system of four spins would require 
15 cyclic permutations, but a method for achieving this with a weighted sum 
of only five permutations has been described |349j . For very large spin sys- 
tems randomly chosen permutations can be used to prepare approximations 
to pseudo-pure states [55] . 

An alternative approach to temporal averaging is to decompose the desired 
pseudo-pure state in terms to product operators. For example the deviation 
density matrix for the pseudo-pure ground state of a two-spin system (Eq. [18]) 
contains the three product operators J^, Sz and Each of these can be 

prepared as the starting state for an individual experiment, and the final result 
obtained by averaging over these starting states as before. Clearly there are 
2" — 1 product operators needed to assemble a pseudo-pure state in an n 
spin system (the spin system is described by 2" product operators, but is 
not required), suggesting that 2" — 1 separate experiments will be required, 
although in practice this can be reduced by combining experiments (obviously 
Iz and Sz can be "prepared" simultaneously). In some cases it can be shown 
by pre-calculation that it is not necessary to include the contributions from 
certain initial states, as these cannot lead to any detectable NMR signals |105] , 
but clearly this approach cannot be used in general. 



4-5 Spatial averaging 

Spatial averaging preceded temporal averaging as a method for preparing 
pseudo-pure states, but has generally proved less popular, perhaps because it 
requires magnetic field gradients which are not always available. Conceptually 
it can be related to temporal averaging much as gradient coherence pathway 
selection is related to phase cycling, although this similarity is not always 
obvious in particular implementations. The original approach of Cory et al. 
[7f8] uses crush gradients to remove undesirable off-diagonal terms, leaving the 
desired diagonal density matrix. It is most easily understood using product 
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operators 
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(53) 




couple 



crush 1 



■>2/z + 25*2 + \2IzSz 



where "couple" indicates evolution under the spin-spin coupling Hamiltonian 
nJ2IzSz for a time 1/2 J. The coupling period can be implemented by refocus- 
ing the Zeeman interactions with spin echoes, but it is not in fact necessary 
to do so, as the initial state commutes with Sz and evolution under the Iz 
term can simply be tracked and then implemented by adjusting the phase of 
the subsequent pulse. In general the action of the gradient pulses in crushing 
off-diagonal terms (which are sensitive to the phase of the reference frame) 
can significantly simplify the implementation of the transformations needed 
to prepare desired states. 

As with any use of gradient crush sequences it is necessary to guard against in- 
advertent gradient echoes, where two crush sequences partly cancel each other 
causing undesired terms to be refocused. In homonuclear systems it is also 
important to ensure that the off-diagonal terms do not include zero-quantum 
coherences, as these are not removed by crush gradients. In heteronuclear 
systems this is not a problem, and the simpler sequence |350] 



can be used instead. It is, however, necessary to begin by equalising the po- 
larisations of the I and S spins; this can be achieved either by applying a 
pulse with an appropriately chosen flip angle to the spin with higher polariza- 
tion followed by a crush gradient, or by using a more complex sequence |350] 
which averages the two polarizations. Sequences have also been developed for 
systems with larger numbers of spins |351j . but these have not been widely 
used. 

Spatial averaging is in some ways preferable to temporal averaging in that 
the result is obtained in a single experiment, but suffers from the obvious 
disadvantage that the signal intensity will be lower. This is not only because 
temporal averaging combines data from multiple experiments, but also be- 
cause the permutation method automatically produces the pseudo-pure state 



Iz + Sz 
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with the highest possible purity, while spatial averaging methods frequently 
sacrifice some of the available purity for simplicity in the preparation sequence. 
One exception to this is use of controUed-transfer gates [352j . which permit 
the spatial averaging equivalent of permutations; related ideas have been ex- 
plored using line-selective pulses |353] . It is also possible to combine ideas from 
temporal and spatial averaging |354] . 



4.6 Logical labelling 



Logical labelling p^fTT] is the conceptually most elegant approach to generat- 

ing pseudo-pure states, although it is quite complex in practice and has not yet 
found widespread use. The basic idea is to find a subset of spin states with a 
pattern of populations matching that of a pseudo-pure state, and then simply 
"relabel" these. For example, a pseudo-pure state of a two qubit system has 
one state with a large population and three other states with equal smaller 
populations, and these can be mapped onto the ground state |000) and the 
three excited states |001), |010) and |100) of a homonuclear three-spin system. 
More sensibly one should choose the ground state |000) and the three upper 
excited states |110), |101) and |011), as this gives a higher signal intensity in 
the pseudo-pure state. Clearly one must always embed the logical spin system 
within a larger number of physical spins, but the required overhead is not 
large [TO] . 

As described the preparation sequence is extremely simple, comprising nothing 
more than mental relabelling of states, but this apparently simple approach 
would come at a great cost in the complexity of implementing quantum logic 
gates, as gates which act in a simple way on the logical qubits will be very 
complex when applied to the physical qubits. A far simpler and more sensible 
approach is to permute the initial populations, so that the desired population 
pattern is moved to the four states |000), |001), |010) and |011), giving a direct 
relationship between logical and physical states pUfTT] . In this case the spin- 
system is in a pseudo-pure state, conditional on the first qubit being in state 
|0), and logic gates can now be implemented directly as long as they do not 
interchange states in the |0) and |1) manifolds of this labelling spin. Explicit 
methods for doing this have been described [11]. A key disadvantage of logical 
labelling remains, however, in that it is necessary to "waste" at least one spin 
as a labelling spin. Applications to date have concentrated on encoding two 
logical qubits in a three-spin system |^95|355] . 




Finally there has been interest in combining logical labelling with averaging 
techniques, allowing the simple preparation of a pseudo-pure state of n — 1 
qubits in an n spin system. The labelled temporal averaging process is usually 
described as labelled flip and swap [SS], and the spatially averaged logical 
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labelling technique (SALLT) is broadly similar |356] . 



^.7 Other routes to pseudo-pure states 

In addition to these three main approaches, many other techniques have been 
described for preparing pseudo-pure states, or approximations to them. Al- 
though superficially different from one another, many of these can be related 
to each other and to the approaches described previously |357] . One interesting 
approach is based on selecting highly entangled states sometimes known 
as cat states after Schrodinger's cat. Such states can be easily generated from 
pure states using simple networks, as illustrated 



for a system with four qubits. Such cat states are closely related to, although 
not quite the same as [T9|2T] . maximal quantum coherences in conventional 
NMR; in particular they have the same n-fold sensitivity to frequency offsets 
and to phase shifts, and so gradients or phase cycling can be used to select 
them from complicated mixtures of states. (This n-fold sensitivity of cat states 
is also the basis of entanglement assisted magnetic field sensing |114II115] .) 

If the cat state preparation sequence is applied to a highly mixed state, such as 
a thermal equilibrium state, then the component in the desired ground state 
will be converted to a cat state, while other components will be converted 
to states with different coherence-orders. Selecting for the cat state is then 
equivalent to selecting the ground state component, and reversing the cat state 
network will generate the desired ground state; the remaining components are 
converted into maximally mixed states, and so this process prepares a pseudo- 
pure ground state. 

Unfortunately it is not possible to apply this procedure exactly as described: 
although cat states do uniquely show an n-fold sensitivity to phase shifts, 
the sensitivity is the same for the two cat states (|0) ± \l))/y/2. Selecting 
for both these states is equivalent to selecting both the desired ground state 
and a second population state, corresponding to a pseudo-pure state of n — 
1 spins in an ra-spin system with the remaining spin being in a state 

proportional to a^- In some cases this apparent waste of a single spin is not, 
in fact, problematic [71] as the mixed state spin can still be used. Despite 
this disadvantage, generation of maximum quantum coherences provides a 
natural route to the preparation of pseudo-pure states, and this approach 
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has been explored in the context of dipolar coupled systems |132fl36] . It is in 
general much simpler to produce states which are almost pseudo-pure than full 
pseudo-pure states and several techniques for doing this have been described 
|13Uf358f359] . These ideas have also been explored in liquid crystal systems 
fT32]|133yi34|[T35]. 



Instead of preparing pseudo-pure states, it is possible to select a pseudo-pure 
component during the final detection of the NMR signal. In some cases the 
signal from the desired pseudo-pure state is localised to a single line in a 
multiplet |114|115] . making its selection particularly simple. A more complex 
approach is to use state tomography, described in Section 15. 3^ to completely 
characterise a mixed state, and then use an eigenvector analysis to select the 
component corresponding to the ground state [76]. 



4-.8 Non-Boltzmann states 



The approaches described above are all methods for preparing pseudo-pure 
states from a thermal mixed state in the high temperature limit [kT ^ hu), 
as this is only regime accessible at equilibrium in the liquid state. To move 
beyond this it is necessary to prepare non-equilibrium states, in which the spin 
temperature is very different from the bulk temperature. Here I will consider 
techniques for generating low temperature spin states, so that kT < hv. Note 
that unless the spin temperature is very low, it will normally still be necessary 
to prepare a pseudo-pure state, but high polarization pseudo-pure states can 
in principle sustain entanglement and so such methods can be scalable. 

Although there are a wide range of conventional NMR techniques for enhanc- 
ing spin polarisation |20|360|361|362|363j . very few of these are capable of 
reaching the low temperature regime in the liquid state. The two main candi- 
dates for doing this are optical pumping and the use of para-hydrogen; a third 
possibility is the use of algorithmic cooling to cool a small subset of spins in 
a larger spin system. 



4-9 Optical pumping 



Spin-exchange optical pumping of noble gas nuclei |364] . notably ^He, ^^^Xe 
and ^^^Xe has been used in a variety of NMR experiments, such as NMR 
imaging |365] . The process involves optical pumping of the electronic states of 
alkali metal atoms, followed by spin exchange between the electron and nuclear 
spin states. Noble gases are, however, unsuitable for NMR quantum computers 
as they are isolated atoms, and so can only provide single qubit systems. It is 
in principle possible to transfer the polarization to more complex molecules by 
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cross polarization |366] , but although this approach has been used to enhance 
signal intensities in NMR quantum computing |367] the transfer efficiency is 
too low to be genuinely useful. Optical pumping in bulk semiconductors |368] 
is even less likely to be useful for liquid state NMR computers. 



4.10 Para-hydrogen 

Para-Hydrogen Induced Polarization (PHIP) |336II337II338] is the only tech- 
nique which has so far been used to prepare NMR spin systems with purities 
above the entanglement threshold [50j. The properties of para-hydrogen are 
a consequence of the Pauli principle, which requires that H2 molecules in 
even rotational states have an antisymmetric nuclear wave function, that is 
the unique nuclear spin singlet state, while molecules in odd rotational states 
have nuclear spin triplet states. Preparation of pure para-hydrogen is diffi- 
cult because interconversion of singlet and triplet states is forbidden; however 
adsorption onto a catalyst surface breaks the symmetry, allowing intercon- 
version to occur. Upon moving away from the surface interconversion is once 
again suppressed, and so the molecular spin state effectively remembers the 
temperature of the last catalytic surface encountered, and remains stable at 
room temperature. A preparation temperature around 20 K is low enough to 
prepare hydrogen with an essentially pure nuclear spin state [50] . 

The para-hydrogen molecule cannot be observed directly by NMR because of 
its high symmetry, but it can be made to undergo a chemical reaction, pro- 
ducing a new molecule in which the two hydrogen atoms are distinct and can 
be observed and addressed. In conventional PHIP experiments this reaction is 
fairly slow, taking around 1 s, and the off-diagonal terms in the initial density 
matrix dephase at a rate determined by the frequency difference between the 
two spins, thus converting a pure nuclear spin singlet state into a mixture 
of singlet and Tq triplet states. In order to prepare a pure ground state for 
NMR quantum computing it is essential to suppress this dephasing. One pos- 
sibility is to apply an isotropic mixing sequence during the reaction 09], but 
relaxation during the reaction remains a problem. 

The most direct approach is to ensure that the reaction is rapid compared with 
dephasing and relaxation, and this can be achieved by using a rapid addition 
reaction with an unstable intermediate prepared by flash photolysis |SU] . This 
leads to an essentially pure spin singlet state which can easily be converted to 
a pure ground state |00), and used for quantum computation [50|51f52f369] . 
So far, however, the process has only be used to prepare two-qubit computers, 
although this polarization could in principle be shared among three qubits 
[370] while remaining above the entanglement threshold. 
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4-11 Algorithmic cooling 



Algorithmic cooling refers to a family of techniques which use computational 
ideas to cool part of a spin system, transferring the undesired heat either to 
other parts of the spin system, or to the wider surroundings. The original idea 
|371j is based on the observation that a molecular heat engine can be made 
to operate reversibly, and so can be implemented using only unitary quantum 
logic gates. Similar ideas were previously considered in the context of polarisa- 
tion transfer |361j . but by using the full control offered by universal quantum 
logic gates it is possible to surpass the S0rensen bounds, and generate spins 
with arbitrarily high polarisation. The key step in this process, a fredkin 
gate, or controlled-SWAP, has been demonstrated in a three-spin system [372j. 

Although these methods do in principle allow arbitrarily high polarisations 
to be reached, the overhead needed is extraordinarily high |371] . with about 
10^ thermal spins being needed to produce a single fully polarised spin |372j . 
Fortunately more practical approaches have been developed, with smaller over- 
heads. Algorithmic cooling |373j surpasses Shannon's spin entropy bound by 
permitting the hot spins to rethermalise with a lattice, so that the entropy 
is now conserved over the whole system and not just the spins in the quan- 
tum computer. To work efficiently this process requires that the spins at the 
hot end of the molecular heat engine have much shorter Ti times that the 
spins which are being cooled. The basic ideas have been demonstrated in ex- 
periments [120|122j . and more sophisticated algorithms have been developed 
|374|375|376|377|378j . but to date genuinely useful polarization enhancements 
remain out of reach. 



4-12 Use of mixed states 

The simplest approach to the problem of preparing pseudo-pure states is just 
to sidestep it completely by using mixed states directly. One example is entan- 
glement assisted magnetic field sensing, where it is not necessary to begin by 
preparing a pseudo-pure state as the desired component can be immediately 
identified as the signal from one of the outermost lines in a multiplet and the 
remaining lines can also be interpreted in a useful way |114fll5] . 

A more generally useful approach is to note that in many computations it is 
not necessary to prepare all the qubits in pseudo-pure states; instead some 
qubits can be left in mixed states [379]. Indeed it has been shown |380] that 
useful quantum computations can be performed in systems with a single pure 
qubit, with the other qubits being in maximally mixed states, an approach 
known as DQCl (deterministic quantum computing with one quantum bit). 



46 



Although DQCl is less powerful than a general purpose quantum computer, 
it can perform some tasks exponentially faster than the best known classical 
algorithm |381f382f383f384] . If the DQCl model is expanded to permit pro- 
jective quantum measurements on a single qubit then it can even be used to 
implement Shor's quantum factoring algorithm |385] : however this is less im- 
pressive than it might seem as the ability to perform projective measurements 
on a single qubit, when combined with quantum swap operations, permits a 
pure state of the whole system to be prepared. 



5 Readout 



The final stage in a quantum computation is reading out the answer, which 
requires the final quantum state to be characterised. This would normally be 
achieved by performing quantum measurements, which project the quantum 
state onto the measurement basis. This approach is not possible in conven- 
tional NMR experiments, as detection of the NMR signal is not a strong 
projective measurement, but rather is a weak ensemble measurement |386] 
which monitors the state of the spin system without changing it. This is easily 
seen by noting that observing the free induction decay in an NMR experiment 
is equivalent to continuously monitoring two non-commuting observables, ax 
and ay. The explanation for this behaviour lies in properties of the operator 
corresponding to the mean of some observable evaluated over an ensemble of 
N identical subsystem; in particular it can be shown that in the limit N oo 
every state is an eigenstate of every mean operator [386J. For finite values of 
N the accuracy of this approximation grows with a/ZV, and for NMR systems 
the error is negligible. 

This ability to monitor without disturbance might seem to be an advantage 
of NMR over other, non-ensemble, approaches, but is in fact a source of many 
problems. In particular, measuring a qubit is the best way of initialising it 
as discussed in Section HI It is, however possible to demonstrate many quan- 
tum algorithms and phenomena using only weak measurements. Just as for 
initialisation, a range of approaches have been used, but all these approaches 
combine the same main elements, based around the analysis of NMR spectra. 
For simplicity I will start by assuming that the computation ends with all the 
qubits in eigenstates, rather than in superposition states or entangled states; a 
small number of measurements will then provide all the information required. 
A more thorough approach is to characterise completely the final state of the 
spin system by quantum state tomography; while this can be useful in small 
spin systems the effort required for full tomography increases greatly with the 
size of the spin system 
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5. 1 Analysis of spectra 



The simplest situation is a one-qubit NMR quantum computer for which the 
states |0) and |1) correspond to the NMR states and — respectively. A 
90° Jy pulse will convert these to ±/a,, which will appear as absorption or 
emission lines in a properly phased NMR spectrum; this makes sense as |0) 
and |1) correspond to excess population in the low energy and high energy spin 
states. As the absolute phase of an NMR signal is meaningless it is essential to 
obtain a reference signal against which phases can be determined; this is most 
easily achieved by acquiring a signal from the computer when it is known to 
be in the initial state |0). 

The situation is similar with larger spin systems, but it is necessary to be 
more careful. The NMR spin state corresponding to |00) is not + Sz but 
h+Sz+'^hSz- A general pseudo-pure eigenstate of two qubits can be expressed 
similarly as 

\ab){ab\ = i ((-!)% + (-l)''^. + 2/,5,) . (56) 

This state can be analysed in two ways. In a homonuclcar spin system it is 
simple to excite and observe both spins, and the two population terms (/^ 
and Sz) are converted to observable single quantum coherences, while the 
longitudinal two-spin order {21 ^8^) becomes unobservable double and zero 
quantum coherence. The observable signal after a 90° pulse is then 

{-lYl, + {-lfS, (57) 

and the desired state information is encoded in the phases (absorption or 
emission) of the NMR signals from the two spins. In a heteronuclear spin 
system it is more natural to excite and observe just one spin, say /. Application 
of a 90° ly pulse to the state gives 

((-1)% + (-1)"^, + (-1)"+^ 2I,Sz) /2, (58) 

and the observable signal is now proportional to 

{-lf{h + {-lf2hSz). (59) 

Only one of the two lines in the / spin doublet will be observed; the choice of 
line depends on 6, the state of spin 5, while the phase of the signal indicates 
a, the state of spin /, as before. 
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5.2 Forced decoherence 

The effect of a projective measurement on a qubit is to project a superposition 
of tfie form 



onto its eigenstates, so that qubit is found in |0) with probabihty |a;p and in 
|1) with probabihty If the outcome of the measurement is unknown then 
the final state is described by a mixed state of the form 



Thus the key effect of projective measurement is to decohere the state, remov- 
ing the off-diagonal coherence terms [387]. 

The same result can be achieved without explicit measurements by forcibly 
increasing the decoherence rate. This can be conveniently simulated by re- 
ducing the apparent coherence time T| by applying a magnetic field gradient. 
Of course gradients do not truly decohere the state, as their effects can be 
reversed by echoes, but diffusion within the sample volume means that gradi- 
ents cannot be exactly reversed |388] . and so gradient dephasing can be made 
indistinguishable from true decoherence [77] . 

Forced decoherence of this kind is widely used when it is desirable to sim- 
ulate the effects of true projective measurements. It can also be useful to 
remove undesired off-diagonal terms which arise from pulse sequence errors 
[44] . although in homonuclear systems it is necessary to worry about the in- 
vulnerability of zero-quantum error terms to gradient dephasing [H]. This 
approach can be generalised to select other states; for example singlet states 
can be selected by twirl operations [229] . Methods for simulating arbitrary 
decoherence processes have been described [5S9j . 

Artificial decoherence can also be imposed using temporal averaging rather 
than spatial averaging approaches, by averaging the result over experiments 
in which a range of gates are applied to the system [8^102] . This approach 
allows considerable control over the exact form of the artificial decoherence, 
and also permits the decoherence to be applied rapidly within a particular 
instance of the experiment, although the use of temporal averaging means 
that the whole process takes longer than spatially averaged approaches. 

5. 3 Tomography 

Quantum state tomography refers to the process of completely characterising 
the state of a quantum system. Originally developed in optics [390] , the process 
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is comparatively simple in spin systems where it is only necessary to determine 
a finite (if large) number of elements in the density matrix description. Clearly 
this is not possible in a single experiment, and tomography can be quite a 
lengthy process. 

In an NMR context it is most helpful to consider the process in terms of 
product operators |116j . A single spin is described by four operators, ^E, I^, ly 
and Iz- Two of these, and Jy, are NMR observables, and can be characterised 
directly by recording the free induction decay. The term is not observable, 
but can be characterised by applying a 90° pulse to convert it into one of the 
observable terms. Note that this process inevitably renders one of the directly 
observable terms unobservable, and it is impossible to characterise these three 
terms in a single experiment. The term is completely unobservable in any 
NMR experiment, and these methods only determine the deviation density 
matrix [TT] . It is also necessary to obtain a reference spectrum from a spin in a 
known state (usually thermal equilibrium), firstly to obtain a phase reference 
allowing and ly to be separated, and secondly to obtain an amplitude 
reference so that the absolute size of the terms in the deviation density matrix 
can be determined. 

In a two-spin system fifteen product operators have to be characterised [391 ], 
of which eight are directly observable, the four inphase coherences (J^, ly, 
Sx, Sy) and the four antiphase coherences {21^3^, 2IySz, 2IzSx, 2IzSy). In a 
homonuclear system these eight terms can all be determined from the am- 
plitudes and phases of the four lines in the spectrum. The remaining seven 
indirectly observable terms must be converted to observables, and it is not 
possible to achieve this for all seven terms in one go. A minimum of four 
experiments is necessary |391|392j . although it is common to use a larger 
number. One popular approach is to use nine experiments corresponding to 
three possible pulses (no pulse, 90° and 90°) applied to spin I and spin S in 
all possible combinations |391] . Quantum state tomography has been widely 
demonstrated in spin-i NMR experiments |llfl3ll53f56f64ll96f393] . as well as 
with quadrupolar nuclei |144fl45f394] and electron spins |164] . 

It is possible to go beyond quantum state tomography to quantum process to- 
mography |395f396ll397j . which attempts to characterise the detailed behaviour 
of a pulse sequence by effectively determining the corresponding superopera- 
tor matrix. This process is inevitably extremely time consuming, but has been 
demonstrated in liquid state |59|226j and solid state |153j NMR experiments. 
Because of the extreme effort required for full process tomography there has 
been interest in partial tomography, which either seeks to determine the over- 
all fidelity of the process [398 1399f400] or concentrates on transfer between 
states of particular spins ^74,.75] . As relaxation is a major source of errors in 
quantum processes there is also considerable overlap between quantum process 
tomography and studies of decoherence. 
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5-4 Methods for single spin detection 



The prospects for NMR quantum computation would be greatly improved if it 
was in fact possible to perform conventional projective measurements on sin- 
gle spins. Direct detection of NMR transitions of a single spin is impractical 
due to the very low energies involved, and it is necessary to couple the NMR 
transition to a process with a larger energy scale. One possible approach is op- 
tically detected ENDOR |401f4U2j . but most interest has focussed on Magnetic 
Resonance Force Microscopy (MRFM) |403f4U4|4U5] , and there have been pro- 




posals for building quantum computers with this approach |406f407f408j . It 

should be noted, however, that such devices will necessarily bear little resem- 
blance to current liquid state NMR experiments. It has been suggested |409] 
that spin-sensing chemical reactions can be considered as projective quantum 
measurements, but this interpretation is controversial. 

A quite different perspective is provided by entanglement based spin measure- 
ments, which seek to enhance the signal from a single spin by correlating it 
with a large number of other spins |41U] . The basic method uses cat state 
techniques to entangle the single spin with a large number of ancilla spins 



and then performs a collective measurement on all the spins. The cat-like state 
corresponds to a state of maximal quantum coherence, and is not directly de- 
tectable, but if the off-diagonal terms in the density matrix are dephased then 
the remaining diagonal terms correspond to population states with very high 
longitudinal spin order, which will give a correspondingly large collective sig- 
nal. This approach has been demonstrated in small systems [228] and possible 
implementations in much larger systems using techniques adapted from quan- 
tum cellular automata have been analysed |332f411f412f413f414] . 



6 Decoherence 

Decoherence |387j refers to any incoherent processes which cause the spin 
system to undergo non-unitary evolution. The term is essentially equivalent 
to relaxation. Some authors chose to distinguish between T2 processes, which 
they call decoherence, and Ti processes, which they call relaxation, but such 
distinctions should not in general be relied upon. As quantum computation 
rehes on the use of coherent superpositions of quantum states, decoherence 
appears to provide a fundamental limit on the complexity of computations 
which can be performed |415] . The development of decoherence free subspaces 
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tolerant computation |425f426f427j have enabled these hmits to be partially 
sidestepped, but decoherence remains a fundamental concern. 



6.1 Modelling decoherence 



Although decoherence in large spin systems can be a very complex process, 
some progress can be made in simulations by incorporating very simple models 
of decoherence [113] , in essence assuming that every spin relaxes independently 
with its own time constants Ti and T2. Typically relaxation is dominated by 
T2 processes, but these largely lead to signal loss, without introducing sig- 
nificant error signals; by contrast Ti processes, although usually significantly 
slower, can lead to ambiguous results in some quantum computations. 

Both types of relaxation are easily modelled using the operator sum represen- 
tation [23] 

p^Y^EkpEl (63) 

k 

where the operators Ek can take any form (in particular, they need not be 
unitary), subject only to the requirement that 

T.ElEk = l (64) 

k 

if the process is to be trace preserving (so that density matrices remain density 
matrices). T2 processes are easily described using phase damping [113] . which 
relies on the fact that for a single qubit the process 

p^lp+^ZpZ (65) 

will completely remove any off-diagonal terms from a density matrix p. More 
realistic decoherence can be modelled by removing the off-diagonal terms from 
some steadily increasing fraction of the density matrix. This gives the set of 
operators 

Eo = ^/x(^^\ Ei = VT^ r " I (66) 





where 

A = l(l+e-*/^^). (67) 

Ti processes are a bit more complicated, as it is necessary to allow for the fact 
that the system does not relax to the ground state, but rather to the thermal 
equilibrium state. This process, generalized amplitude damping [113] . can be 
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described by the set of four operators 
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where 



7 = 1- e-*/^i 



P=2 




(69) 



and e is the polarization at thermal equilibrium, so that the thermal state is 



As expected pth is unaffected by either kind of relaxation. It might seem tempt- 
ing to concentrate on the pseudo-pure component of the spin state, which 
corresponds to e = 1, in which case the operators E2 and can be dropped, 
but while this process will work for a single spin it will give incorrect results 
when extended to larger spin systems. 

The complete relaxation process for a single spin can be described by applying 
both amplitude damping and phase damping, remembering that the measured 
T2 time includes a contribution from Ti relaxation. The relaxation of a group 
of two or more spins can be described by applying relaxation sequentially to 
each spin; in this case it is vital to calculate the amplitude damping correctly, 
so that the system relaxes back to the correct thermal state and not to a 
pseudo-pure state. 

The approach outlined above only applies to decoherence during a single pe- 
riod of free precession. With a more complicated pulse sequence the simplest 
approach is to break it up into a sequence of pulses (assumed to be instan- 
taneous) and delays, with the decoherence that occurs during a delay being 
calculated at the end of the delay, immediately before the following pulse. Note 
that only the spins experiencing the pulse need be considered at this stage; 
decoherence of other spins can be saved up, to be applied just before the next 
pulse addressing those spins. If only T2 processes are considered then this can 
be simplified further, as T2 relaxation is not affected by 180° pulses, and so 
T2 relaxation can be calculated at the end of a whole spin-echo period, rather 
than step by step. This approach will not be entirely correct in experiments 
with long selective pulses, but has proved fairly successful in some cases |113] . 
It will, however, be difficult to apply to pulse sequences which make extensive 
use of gates derived by optimal control methods, as in this case any division 
into pulses and delays is to some extent artificial. 
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6.2 Decoherence in large spin systems 



The previous section has assumed that every spin relaxes individually, uncor- 
related with the relaxation of neighbouring spins, and it is well known that 
this assumption is rarely correct in NMR spin systems. Obvious examples in 
conventional NMR include nuclear Overhauser effects [428] arising from corre- 
lated Ti relaxation, the differential linewidths in multiplets used in transverse 
relaxation-optimized spectroscopy (TROSY) experiments |429] . and the relax- 
ation of multiple quantum coherences |430j . 

NMR techniques can also be used to study decoherence processes, and in par- 
ticular to determine how decoherence rates scale with the number of qubits 
in real physical systems, as well as the use of decoupling techniques to re- 
duce decoherence rates by suppressing undesirable interactions in such systems 
1106010702640431043204331. 



7 Quantum Algorithms 



I will now consider some of the quantum algorithms which have been imple- 
mented using NMR quantum computers. Before doing so it is useful to say 
a little more about the implementation of mathematical operations with re- 
versible quantum devices, and in particular to explore how functions can be 
evaluated in a reversible context. The obvious direct approach of creating a 
propagator which implements 

\x) % \f{x)). (71) 

is not normally practical, as this process will only be reversible if the function 
is itself intrinsically reversible, which is not generally the case. The solution 
to this problem is to preserve the input qubits and store the output in an 
ancilla bit. The output cannot, however, simply overwrite the initial value of 
the ancilla, as that would be irreversible. Instead the new value of the ancilla 
must be obtained by reversibly combining the output of the gate with the old 
value, and this is most simply achieved by using bitwise addition modulo 2 
(the XOR gate), so that 

\x)\y)^\x)\y(Bf{x)). (72) 

If the ancilla bit is initially set to zero then the value of the function can be 
easily obtained 

|x)|0)^|x)|/(a;)). (73) 
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For a function with two input bits and one output bit this can be written as 
a quantum network 

(74) 



b — 



f 



— b 

- c® f{a,b) 



c ^ 

sometimes caUed a /-controUed-NOT gate. Functions with more than one out- 
put bit can be handled by combining each output bit with its own ancilla. 
Note that while reversible logic can perform any desired transformation on 
a set of input bits it does not provide any means to set the input or ancilla 
bits into the desired initial states, and so a reversible computation must begin 
with an irreversible initialisation process. 



7.1 Deutsch's algorithms and related methods 



Deutsch's algorithm plays a central role in the historical development of quan- 
tum computing, being both the first quantum algorithm to be described and 
the first quantum algorithm to be implemented experimentally. It considers 
the analysis of the simplest type of functions: binary functions which map 
a single input bit to a single output bit. Clearly there are only four such 
functions 
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where each function is conveniently labeled by the pattern of output bits in 
its truth table. 



Deutsch's problem considers the identification of a function, which is known 
to be one of these four functions but otherwise undetermined, in a situation 
where the only possible access to the function is by finding its output for a 
particular input (that is, there is no way of determining how the function val- 
ues are calculated). This is sometimes called the oracle model of computation, 
and each function evaluation is called an oracle query. Clearly it is necessary 
and sufficient to make two oracle queries to completely identify a function, 
evaluating both /(O) and /(I) with /-controUed-NOT gates. 

The function can in some cases be partially identified in a single query: for 
example /oo and /oi can be distinguished from /lo and fu by evaluating /(O). 

Suppose, however, it is desired to distinguish functions according to their 
parity, /(O) © /(I), so that the balanced functions /oi and /lo, which give 
outputs of for one input and 1 for the other, are distinguished from the 
constant functions /oo and fu, which give the same output for all inputs. 
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With a classical computer this distinction requires two queries, as calculating 
the parity requires both /(O) and /(I) to be determined. Deutsch, however, 
showed that a quantum computer can determine the parity with a single query. 

Deutsch's original algorithm is not particularly useful, as it only works half 
the time at random [1], reporting an inconclusive result the other half of the 
time. Thus references to Deutsch's algorithm almost always refer to a more 
modern version which works with certainty [434] . The quantum network to 
achieve this 
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is well explained elsewhere p^21|l24|l434] . but can be broken down into three 
parts. The initial Hadamard gate on the first qubit generates |+), a uniform 
superposition of the two possible input values; this is then passed to the /- 
controUed-NOT which evaluates the function for both inputs simultaneously. 
As the ancilla qubit is in the state H|l) = |— ) the effect of the function on a 
single input is given by 
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and this result can be simplified by considering the two possible values of f{x). 
For the case f{x) = the result is 



|x)|0)-|x)|l) 



V2 



1^)1-) 



(7J 



while if f{x) 



1 the result is 
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This can be summarised as 
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so the value of the function determines the phase of the state. This might 
appear to be a global phase, and so irrelevant, but if the input qubit is in a 
superposition state then the function values can appear as a relative phase 
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If the function / is constant, so that /(O) = /(I), this simplifies to 

-) = (-l)^^°^l+)|-) 
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while if the function is balanced the result is 



(_1)/(0)MH1)|_^^(_^)/(0)|_^|_^_ 



^3) 



The effect of the final pair of Hadamard gates is to convert this to |0)|1) or 
as appropriate, completing the desired computation of the parity. Note 
that the final value of /(O) is encoded as an undetectable global phase, and so 
this algorithm enables only the parity to be determined; it is not possible to 
make further distinctions between the two constant or two balanced functions 
without making further oracle queries. 

The simplicity of the network, requiring only two qubits, makes it an extremely 
attractive early target. It is necessary to implement the /-controlled-NOT gates 
explicitly, but these are trivial for the two constant functions, /oo and /n, and 
only require a single controlled-NOT gate for the two balanced functions, /oi 
and fiQ [T^pi] . Deutsch's algorithm was implemented almost simultaneously 
using homonuclear [12] and heteronuclear [13| two-spin systems, and has sub- 
sequently been implemented on an ion trap quantum computer p^5j, with a 
pure state NMR quantum computer [51], and with a photon cluster state 
quantum computer |435] . 

Most NMR developments have, however, involved not Deutsch's algorithm, 
but its generalisation, the Deutsch-Jozsa algorithm |434II436] , which considers 
functions from n bits to 1 bit (as Deutsch's algorithm can be considered as the 
simplest example of the more general Deutsch-Jozsa algorithm, with n = 1, 
the general name is sometimes used for this well [19|63] ) . Such functions 

need not be constant (the same output for all inputs) or balanced (an equal 
number of Os and Is in the outputs); for example a function with n = 2 
might output for one input and 1 for three inputs. Suppose, however, that 
the function is promised to be either constant or balanced; in this case it is 
possible to determine which sort of function it is with a single evaluation. The 
quantum network is an expansion of that for Deutsch's algorithm, and for 
n = 2 takes the form 
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where the output bits, a and b, are both equal to if the function is constant, 
and are not both equal to if the function is balanced. This generalises to n 
bit functions in the obvious way, with all output bits equal to if the function 
is constant, and at least one bit equal to 1 if the function is balanced. 

Implementing the Deutsch-Jozsa algorithm would appear to require a device 
with at least three qubits, and the first implementation p9] did precisely that, 
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using line selective pulses in a three-spin homonuclear spin system; this simpli- 
fies the implementation of two of the balanced functions, which would require 
TOFFOLI gates if implemented using standard circuits |355f48flU4] . The algo- 
rithm has also been implemented for 4-bit functions in a 5-spin system |lU5j . 
restricting the choice to a single balanced function which could be implemented 
relatively simply. 

Although the conventional Deutsch-Jozsa algorithm requires a minimum of 
three qubits, it was noted early on |437] that the ancilla bit is not really 
necessary. The ancilla bit is used in classical function evaluations to hold the 
function value, but in quantum versions the effect of setting the ancilla bit to 
|— ) is that the function value is stored as a phase rather than a bit value. Thus 
it is possible to remove the ancilla qubit if the conventional classical oracle is 
replaced by a purely quantum oracle, which implements 

\x) ^ (-l)^W|a;) (85) 

directly. This quantum oracle cannot be used for classical algorithms (as the 
qubit phase is a purely quantum mechanical phenomenon), but will work 
perfectly well in quantum algorithms. This purely quantum version [SS1S2], 
which is sometimes called the refined Deutsch-Jozsa algorithm [S3] or the 
Collins version of the Deutsch-Jozsa algorithm [96], permits ?7,-bit functions 
to be studied using an ra-spin system. This approach can even be used to 
implement Deutsch's algorithm with a single spin |47j . 



7.2 Grover and other quantum search algorithms 

Grover's quantum search algorithm [43811439] is the simplest member of the 
second major family of quantum algorithms. In general the algorithm concerns 
the analysis of functions of n input bits for which the output is except 
for a single satisfying input, for which the output is 1. If the function is 
only accessible through an oracle implementation, so that it is not possible 
to analyse the function itself, but only its values, the only practical approach 
to locating the satisfying input is to try inputs at random until it is found. 
This unstructured search, equivalent to trying to find the name corresponding 
to a particular telephone number by reading through the telephone directory, 
is obviously extremely inefficient: for a function with n input bits there are 
N = 2^ possible inputs, and on average it will be necessary to try about N/2 
of these before the satisfying input is found. By contrast, Grover's algorithm 
allows it to be located with around a/ZV queries. 

The quantum circuit implementing Grover's search algorithm in the case n = 2 
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{N = 4) is 
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where the first two bits represent the inputs to the function and the third bit 
is an ancilla. The first three qubit gate is an /-controlled-NOT gate and the 
second one is a close relative of the toffoli gate which applies a not gate to 
the ancilla bit if the input bits are both in state 0. As in the Deutsch-Jozsa 
algorithm the ancilla bit is not really necessary as it only serves to convert 
the function result from a bit value to a phase, and so the network can be 
simplified to 
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where Uqq negates the basis state |00) but leaves all other basis states alone, 
and Uf negates the basis state corresponding to the unique satisfying in- 
put. This simplified algorithm is universally used in NMR implementations 
of Grover's search. 



As usual the initial Hadamard gates convert the input qubits to a uniform 
superposition of the four possible inputs, and the function is evaluated over 
all four inputs in parallel, leaving the satisfying input uniquely marked with 
a phase shift. The remaining gates then convert this phase difference to an 
amplitude difference, so that all the amplitude in the superposition is concen- 
trated on the satisfying input, and when the input bits are measured at the 
end of the process their values will indicate the satisfying input. The circuit 
is similar for cases with n > 2, except that it is now necessary to repeat the 
function evaluation and amplitude concentration stages several times, with 
the amplitude of the satisfying input in the final superposition varying sinu- 
soidally with the number of repetitions. This has two consequences: firstly 
that it is necessary to choose the number of repetitions appropriately, and 
secondly that it is usually not certain that the satisfying input will be located, 
even if the optimum number of repetitions are used, as in the general case the 
amplitude may never be completely concentrated on the satisfying input. 

The case n = 2 is special, in that the algorithm is guaranteed to work af- 
ter a single repetition, making this case an attractive target, and Grover's 
search was implemented early on with both heteronuclear [31] and homonu- 
clear [13] spin systems. Later implementations have been used to demonstrate 
the use of dipolar couplings in liquid crystal solvents [56], efficient quantum 
state tomography [61], the stability of Grover's search under imperfect gates 
[67] , and NMR quantum computing with highly polarized [367] and almost 
pure [52] states. Unlike the Deutsch-Jozsa algorithm there has been little in- 
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terest in increasing the value of n, but there has been one demonstration of 
searching over three-bit functions on a three-spin device [79], demonstrating 
the expected oscillatory behaviour as the number of repetitions is varied. 

Grover's algorithm can be extended to more general functions with more than 
one satisfying input; the network is identical, but the number of repetitions 
required now depends on ^J N/k, where k is the number of satisfying inputs 
|440] ■ The output qubits are in a uniform superposition of all the matching in- 
puts, which makes this an unattractive algorithm for NMR implementations, 
as the non-pro jective readout process will not cope well with this superpo- 
sition. A related algorithm, approximate quantum counting, permits k to be 
estimated by observing the oscillation frequency |44Uf441] , and this is a much 
more attractive target. The algorithm requires an additional control qubit in 
addition to the n bits defining the input to the function, but a two-spin sys- 
tem has been used |44j to demonstrate the case n = 1, for which k equals 0, 
1 or 2, and a related experiment has been performed in a three-spin system 
[93]. The oscillatory behaviour inherent in quantum counting allows this algo- 
rithm to be used to explore the effects of applying very large numbers of logic 
gates, and this has been used to demonstrate the effects of composite pulses 
in suppressing off- resonance errors [211] . 

An alternative approach developed by Grover replaces the simple repetition 
of a single search operator with a more complicated recursive operator [442] , 
which causes the algorithm to converge on a superposition of satisfying inputs 
no matter what the value of k is. This algorithm permits an efficient search to 
be performed when the value of k is unknown, although it is much less efficient 
than the original quantum search algorithm when k is known. This algorithm 
has been implemented by NMR [61J for the case n = 2 in a two-spin system. 
Results were demonstrated for functions with k = 1 and = 2; in the latter 
case care was needed in interpreting the final result, as this corresponds to a 
superposition of the two satisfying states. 

In addition to Grover's unstructured search, quantum algorithms have also 
been developed for searching functions whose structure is partly known. NMR 
implementations include the Bernstein- Vazirani algorithm [75|443] and Hogg's 
structured search [73j444ll445j . Several of these implementations build on the 
concept of searching in Liouville space [183f446ll44 7] which is particularly ap- 
propriate to ensemble implementations such as NMR. 



7.3 Shor type algorithms 

Shor's quantum factoring algorithm [448f449ll450ll451] is perhaps the best known 
application of quantum computation because of its potential significance to 
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wider society. Many current cryptographic schemes underlying the security of 
financial systems |452j and of the internet in general |453] rely on the appar- 
ent difficulty of factoring large composite numbers, or closely related problems 
such as calculating discrete logarithms |454] . No efficient algorithm is known 
for solving these problems on a classical computer, and it is widely presumed 
that no such algorithm exists. By contrast a quantum computer can efficiently 
solve such problems, rendering much of current cryptography obsolete. 

Shor's algorithm works by finding the period of the modular exponentiation 
function, after which results from classical number theory allow a variety of 
particular problems to be solved |449] . Period finding is, in fact, a general 
feature of several quantum algorithms, such as the Deutsch-Jozsa algorithm, 
and Shor's algorithm can be related to these [450] . A key element is the imple- 
mentation of a quantum Fourier transform [449j which was demonstrated by 
NMR in a three-spin system [85] and used to find the order of a permutation 
in a five spin system |109] . 

It was then a small step to a full implementation of Shor's algorithm by NMR 
|113] . A key step in achieving this |455] was that realisation that the imple- 
menting the simplest non trivial example of the algorithm, in effect factoring 
= 15, is particularly simple, requiring only seven qubits. The implemen- 
tation was further simplified by noting that certain gates could be, in effect, 
compiled out of the quantum network, substantially reducing the complexity 
of the final circuit. Even so the final pulse sequence |113] . containing around 
300 distinct shaped pulses, remains one of the most complex NMR experiments 
ever demonstrated. 



7.^ Quantum simulation 

Quantum simulation is the idea of using one quantum mechanical system to 
simulate another quantum mechanical system more efficiently than can be 
achieved by any system based on classical physics J5j. Although such simu- 
lations may not look like conventional algorithms, a quantum simulator can 
equally well be thought of as a special purpose quantum computer for solving 
an otherwise intractable problem. 

It is wise to use some caution in developing this argument, as in principle al- 
most any physical process can be considered as a computation |456|457] . The 
approach is, however, clearly useful as long as it is possible to simulate evo- 
lution under arbitrary Hamiltonians, and it has been shown |458f459f460f461] 
that any quantum system evolving under a local Hamiltonian can be efficiently 
simulated on a quantum computer. Specific algorithms have been developed 
for many important problems, including molecular energies |462] . molecular 
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geometries |463] , and chemical dynamics |464j . Arbitrary Hamiltonians can ei- 
ther be simulated using quantum networks, or by constructing effective Hamil- 
tonians using Trotter approximations [2^465] . 

Simple simulations can be implemented on a two qubit system, and NMR was 
used early on to simulate a truncated quantum harmonic oscillator [1TP2] . 
Subsequently NMR techniques were used to simulate the Heisenberg Hamil- 
tonian in spin chains |466] . and pairing Hamiltonians |467f468f469] . Recently 
quantum chemistry algorithms developed for use in photonic systems |470j 
have also been applied in NMR |471] . 

7.5 Adiabatic algorithms 

If a quantum system starts in the ground state of its Hamiltonian, and the 
Hamiltonian is slowly varied, then the system will stay in its instantaneous 
ground state. This adiabatic following forms the basis of adiabatic algorithms 
for quantum computing |472|473] . These algorithms work by constructing a 
Hamiltonian corresponding to the problem to be solved, such that the ground 
state of the Hamiltonian encodes the solution. This ground state can then 
be determined by slowly interpolating the experimental Hamiltonian between 
some initial form, such as the background Hamiltonian of the system, and 
the desired final form. Although this approach appears quite different to con- 
ventional circuit methods, it has been shown that the two schemes are fun- 
damentally equivalent |474] . and the methods used to implement adiabatic 
algorithms are closely related to those used in quantum simulations. 

This approach was first implemented using a three-spin system |475] to im- 
plement an optimization algorithm, with the interpolation between the two 
Hamiltonians implemented using Trotter approximations |465] . Subsequently 
adiabatic methods were used to demonstrate Grover's quantum search and 
Deutsch's algorithm ^475], to solve the 1-SAT satisfiability problem pS], and 
to implement quantum factoring |80] . 

7.6 Gauss sums 

Several approaches have been described |477|478|479] for implementing Gauss 
sums by NMR, and it has been suggested that this can be used as a method 
for factoring numbers. This relies on the fact that the Gauss sum 
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(where is an integer to be factored, and / is a trial factor)is equal to 1 if / is 
indeed a factor of N, and has an absolute value less than one otherwise, with 
the exact value depending on M. Several different methods for evaluating this 
sum have been described, but in essence they all involve applying a sequence 
of many pulses to a single isolated spin, with the phase of each pulse effectively 
encoding a single term from the Gauss sum. Related experiments have also 
been implemented with other techniques |480|481] . 

These NMR methods have generated considerable interest, as they can be ap- 
plied to far larger numbers than current implementations of Shor's quantum 
factoring algorithm. It is, however, important to note that Gauss sums do not 
in themselves provide a factoring algorithm, but rather a factor checking al- 
gorithm, and it may be necessary to try all possible prime factors up to \/N. 
Since the implementations only involve a single spin they cannot, of course, in- 
volve entangled states, raising further doubts as to their importance. However 
it has been previously noted that even single qubits using superpositions to 
integrate phase can outperform simple classical systems [482j. Unfortunately 
current NMR methods to evaluate Gauss sums require the initial computation 
of pulse phases, and this cannot be implemented without first implicitly de- 
termining whether or not the trial factor is indeed a factor |483j , and so these 
implementations are not interesting from a computational point of view. It is 
not yet clear to what extent these criticisms also apply to other techniques for 
implementing Gauss sums [484] . 



8 Other quantum information phenomena 

In this section I consider a range of phenomena and techniques from quantum 
information processing which have been demonstrated or explored using NMR 
techniques. As was the case for quantum computation NMR is a useful tech- 
nique for in-principle demonstrations of many of these topics even though the 
lack of projective measurements and pure states means that the techniques 
cannot normally be used for anything genuinely useful. 

8. 1 Quantum communication protocols 

Quantum communication protocols use quantum mechanical effects to trans- 
mit information between two or more parties in a way which is impossible 
with purely classical devices. Perhaps the most important example is quan- 
tum cryptography [485] . which permits two parties, usually called Alice and 
Bob, to securely establish a random key which can be used to encrypt a sub- 
sequent message, but many other examples are known. NMR is peculiarly 
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ill-suited to such applications, not only due to the lack of projective mea- 
surements and pure states, but also because it is effectively impossible to 
spatially separate individual qubits which are encoded as nuclear spins in a 
single molecule, but this has not prevented simple demonstration experiments 
from being performed. 

To begin with, note that a qubit only has two distinguishable states, |0) and 
and so it might seem that only a single bit of classical information could be 
reliably transmitted by a single qubit. (While superposition states appear to 
store more information, in that it takes more than a single bit of information 
to accurately describe them, this information cannot be accessed by measuring 
the state |486j .) Quantum dense coding |487j uses the peculiar properties of 
entangled states to enable two bits of classical information to be transmitted 
by a single quantum bit, which can normally only reliably transmit a single 
bit of information. 

The protocol assumes that Alice and Bob each possess one of a pair of qubits 
in the entangled Bell state and that Alice wishes to transmit a two bit 

message to Bob. She can do this by performing one of four transformations, 
1, X, Y, or Z, to her qubit; as the two qubits are entangled this operation 
affects the whole two qubit state, converting it to lip'), l^"*"), or 

respectively. If Alice then sends her qubit to Bob he can determine which of 
these four Bell states the two qubits are in, and thus which message Alice 
wished to send. This protocol has been implemented with a two qubit system 
|488j . and a generalisation has been demonstrated in a three qubit system 
[i89] . 

Quantum teleportation |490j also uses entangled states, this time to transmit 
an unknown quantum state. Suppose Alice has a qubit in some unknown state, 
which she wishes to send to Bob. As she does not know the state she cannot 
simply send him a description, and it is furthermore impossible for her to 
obtain these details as it is impossible to accurately measure an unknown 
quantum state |486] . Thus it might seem that the only way in which she can 
transfer her state to Bob is to physically transfer the qubit carrying the state. 

Quantum teleportation, however, enables her to achieve this without physi- 
cally transferring the qubit as long as she shares an entangled state with Bob 
and is able to send him two bits of classical information; her original copy 
is destroyed in the process. This protocol was demonstrated in a three-spin 
system [75] , and was indeed one of the earliest quantum protocols to be demon- 
strated by NMR. This experiment relied on rapid T2 decoherence of two ^^C 
nuclei to simulate the effects of quantum measurement, and emulated Bob's 
final state transformation step by post-processing, thus greatly simplifying the 
implementation. 
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Entanglement swapping |487f491ll492] is a close relative of quantum telepor- 
tation, in which the quantum state to be teleported is itself entangled with 
another qubit. As quantum teleportation transfers a perfect copy of the un- 
known quantum state, it will transfer this entanglement. Thus if Alice shares 
an entangled state with Bob, and another such state with Charlie, she can 
teleport her entanglement such that Bob and Charlie now share an entangled 
state. This approach can be generalised to distribute more complex forms 
of entanglement, sometimes called a quantum telephone exchange |493j . The 
scheme has been demonstrated in a four spin system |1U2] , using artificially im- 
posed decoherence to implement the measurement step. Entanglement transfer 
can also be achieved by replacing the teleportation stage with a pair of quan- 
tum SWAP gates [lOlj . but this is not normally considered a communication 
protocol as it requires coherent interactions between the three parties. 

Remote state preparation |494|495] is another generalisation of quantum tele- 
portation, in which Alice wishes to transfer a known quantum state to Bob. 
As she knows the state she could simply transmit an accurate description, but 
this requires transmitting two real numbers, such as the 6 and (p coordinates of 
the state on the Bloch sphere. She could instead teleport the state, requiring 
only two classical bits to be sent, as long as she shares an entangled state with 
Bob. As she knows the state, however, the process can be slightly simplified, 
and the state can be transmitted using a single classical bit. Furthermore it 
is not strictly necessary for Alice to even prepare a qubit in the known state, 
and the whole process can be performed using only the two qubits making up 
the entangled state; this simplified process has been demonstrated in a two 
qubit system |496] . 

Quantum random walks |497f498] are quite different from their classical coun- 
terparts, as a result of quantum interference effects. Both continuous time 
|499j and discrete time |232] walks have been investigated by NMR. 

Quantum chaos occurs when a quantum system evolves under a Hamiltonian 
in such a way that its final state is hypersensitive to small perturbations in 
the dynamics |500] . A simple example, the quantum Baker's map |500j . has 
been studied using a three-spin system [90] . 

Many quantum communication protocols can be recast using the idea of quan- 
tum games |501f502f503f504] , with the participants in the protocol, and some- 
times nature itself, taking the role of players. There have been two explicit 
demonstrations of a traditional quantum game, the prisoner's dilemma, using 
two-spin |5U5] and three-spin |5U6] systems. 
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8.2 Quantum Cloning 



It is a well known result in quantum information theory that a qubit in an 
unknown state cannot be accurately copied {cloned) |507j . It is, however, pos- 
sible to prepare approximate copies |508j . and specific quantum networks for 
doing this have been described. For example an optimal one-to-two cloning 
network |509j converts a single copy of a unknown quantum state into two 
approximate copies of the form 

il^)(^l + |l^^)(^^l = il^)(V^I + |l/2 (89) 

where is the state orthogonal to \ip). This description can be naturally 
generalised to mixed states, and such mixed state quantum cloning has been 
demonstrated by NMR [71]. This implementation used a three-spin system, in 
which the initial spin is cloned onto two ancilla spins. The initial spin is left 
in a mixed state corresponding to an optimal implementation of a transpose 
operation on the original state, which is closely related to a universal-NOT 
gate [5T0II5TT] . 

There are several other processes closely related to quantum cloning, such 
as cloning with prior partial information |512] , which permits more accurate 
copies to be made if the initial state is limited to a subset of states, such 
as states on the Bloch sphere, Eq. [121 with a fixed co-latitude but unknown 
azimuth angle. There has also been an implementation of asymmetric quantum 
cloning |293j in a two-spin system; this process uses no ancilla spins, and allows 
the fidelity of the cloning to be arbitrarily divided between the two copies. 

8.3 Decoherence, error correction, and fault tolerance 

Decoherence processes, that is unwanted and uncontrolled interactions be- 
tween individual qubits, or between qubits and the environment, will lead 
to errors in the values represented by these quantum states. Although de- 
coherence effects can be minimised by careful control of such interactions, 
such errors cannot be entirely avoided, and must be either suppressed or cor- 
rected |513j . The two principal methods for doing this are the use of de- 
coherence free subspaces |416|417|418ll419|420j and quantum error correction 
|421f422|423|424j . These can be thought of as quantum communication proto- 
cols, allowing a quantum state to be accurately transmitted from the present 
to some later time. A further generalisation of this is fault tolerant compu- 
tation |425f426f427] . which considers the effects of decoherence during the 
application of logic gates, particularly those logic gates used to implement 
decoherence free subspaces and error correction. All these methods rely on 
replacing a single physical qubit with a logical qubit, encoded with the help of 
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multiple ancilla qubits. 



When considering methods for handling decoherence it is useful to distinguish 
between two major types: phase decoherence, which can be considered as 
the random application of Z gates to the state, and spin flip decoherence, 
which arises from random X gates; more complex decoherence processes can 
be composed as a combination of these fundamental operations [23], and so it 
normally suffices to protect a system against these two types of decoherence. 
(Note that phase decoherence is essentially equivalent to T2 relaxation, while 
spin flip decoherence is not quite the same as Ti relaxation.) The two processes 
are closely related, as a random X gate is equivalent to a random Z gate 
surrounded by Hadamard gates, and so a method which protects against phase 
decoherence can easily be converted to one that protects against spin flips, 
and vice versa. Protection against both sorts of decoherence simultaneously 
can be achieved either directly or by nesting methods that protect against 
the individual errors; the latter approach is often simpler to implement but 
requires a larger overhead in ancilla qubits. 

Decoherence free subspaces use correlations between the decoherence of in- 
dividual qubits to design multi-qubit states which are (ideally) unaffected by 
decoherence. Although it is common to assume that spins relax independently, 
this is rarely the case as discussed in Section 16.21 above. These effects can be 
considered as implicit demonstrations of decoherence free subspaces, a clear 
example being the use of zero-quantum coherences to observe NMR signals in 
the presence of intense local magnetic field inhomogeneities |514] . These simple 
schemes usually only involve protection of the system against transverse relax- 
ation. Using the theory of decoherence free subspaces [416|417|418|419|420] it 
is possible to design systems which are simultaneously protected against cor- 
related relaxation of all kinds by using multiple spins to encode a single qubit. 

An explicit demonstration has been performed [SI] in which a single qubit 
encoded on three ^^C nuclei was shown to be resistant against artificially im- 
posed decoherence (see Section [5l2|) . and this system was subsequently studied 
in detail |515j . A simpler system in which a qubit is encoded in two spins, pro- 
viding resistance only to transverse relaxation, has also been studied |516] , and 
the use of identical spins in CH„ systems has been explored |517] . However, 
decoherence free subspaces are only really useful when multiple qubits are si- 
multaneously protected. To date a four-spin system has been used to protect 
two qubits against artificially applied spin fiip decoherence jlOOj . and it was 
shown that the resulting device could perform a quantum algorithm, while a 
similar system not using decoherence free subspaces could not. The behaviour 
of systems of this kind has been studied in considerable detail |230f260] . 




While correlated errors can be suppressed using decoherence free subspaces, 
uncorrelated (independent) errors must be detected and corrected. In classical 
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information processing this problem is well understood |176j . and the solution 
relies on redundancy, keeping multiple copies of each bit so that disagreements 
between copies can be identified and reconciled. The simplest approach is to 
keep three copies of every bit; should these three copies ever disagree then 
the single disagreeing bit can be reset to the majority value. This approach 
assumes that the probability of two bits being wrong is negligible, which will 
be reasonable if the errors occur independently and infrequently, so that the 
probability of two errors is much smaller than the probability of one. If the 
error rate is too high for this approximation to be used then more sophisticated 
methods, using larger numbers of copies, can be used. Similarly, if the error 
rate is very low then it is possible to use parity checks to correct errors with 
a smaller number of ancilla bits. 



At first sight it is difficult to see how these ideas could be used in quantum 
computation, as they appear to rely on measuring the bits to detect differ- 
ences between copies, and this cannot be done for bits in superposition states 
without causing the superpositions to collapse. It is, however, not actually 
necessary to measure the bits to detect errors: it suffices to detect whether 
two bits are the same or different, without actually determining their states 
and the erroneous copy can then be rotated back into place. 



Quantum error correction has been implemented using a three-spin system 
[77] to protect against phase decoherence. A related demonstration used a two- 
spin system to implement detection (but not correction) of phase errors |60j . 
permitting a detailed characterisation of the process; a detailed analysis of the 
three-spin system was subsequently performed |518] . Protecting a single qubit 
against arbitrary decoherence requires four ancilla qubits |424] . and this has 
been demonstrated using a five-spin system [112j . Quantum error correction 
has also been successfully combined with decoherence free subspaces using a 
four-spin system |519] . 



All these demonstrations have only implemented a single round of error cor- 
rection. In order to preserve a qubit state indefinitely it is necessary that the 
qubit state be corrected repeatedly, and this in turn requires that the ancilla 
qubits used to detect errors be reset back to their initial state after each round. 
As previously noted, the ability to reset qubits during a computation is not 
available in NMR implementations, and so full error correction is currently 
impossible. Error correction has also been implemented in ion traps |520f521] . 
where resetting qubits is theoretically possible; however as yet no complete 
demonstration of multiple round error correction has been performed. 
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8-4 The quantum Zeno effect 



The quantum Zeno effect |522|I523|I524|I525] refers to the fact it is possible 
to suppress the coherent evolution of a quantum system by making frequent 
measurements, which project the quantum system onto its eigenstates, during 
the evolution process. As an example consider a not gate applied to a qubit 
starting in state |0). (In NMR terms this corresponds to a 180° rotation applied 
to an isolated spin starting in the state I^-) During the evolution the state is 
described by 

\il}{t)) = cos(a;t/2)|0) + isin(a;t/2)|l) (90) 

with nutation rate u = vr/tigo. Suppose, however, that measurements are made 
in the {|0), |1)} basis at time intervals r: each measurement will project the 
system onto either state |0), with probability cos^(a;r/2), or state with 
probability sin^(a;r/2). If measurements are made rapidly, such that ut <^ 1, 
the system will almost always be found after the first measurement in the 
initial state |0). Subsequent evolution and measurements will have the same 
effect, with the system being returned to the initial state each time, so frequent 
measurements can effectively suppress the nutation. 

It is not immediately obvious how NMR can be used to demonstrate a quan- 
tum Zeno effect, as the standard NMR measurement process is not a true 
projective measurement. It is, however, possible to simulate the effects of pro- 
jective measurements by forced decoherence as described in Section [5l2] above. 
This approach has been used to demonstrate the quantum Zeno effect using a 
single spin [22] • The effects of measurement strength were also explored in 
a two-spin system [22] • An analogy can be drawn between the Zeno effect and 
the behaviour of frequency selective pulses [62] , and more generally between 
Zeno effects and dynamical decoupling sequences [526] . 



8.5 Multi-spin entanglement 

Although high temperature pseudo-pure spins states cannot be truly entan- 
gled, there has been interest in preparing pseudo-entangled states [21], that 
is pseudo-pure states whose pure part is entangled. This is sometimes seen as 
a test of quantum computing [111] , as a general purpose quantum computer 
should, among other things, be able to produce highly entangled states in an 
n-qubit system. 

In a two qubit system the maximally entangled states are the Bell states, but in 
larger systems the situation is less simple, as different entanglement measures 
can lead to different conclusions. With three qubits the states considered most 
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frequently are the Greenberger-Horne-Zeilinger (GHZ) states |527f528] 

|^GHz) = (|000)±|lll))/V2 (91) 

and the W state 



\i)w) = (|001) + |010) + |100))/V3 



(92) 



although a more complete classification of states is available |53U] . GHZ states 
were produced in early experiments [76j, although in this case a pseudo-pure 
state was not initially prepared, with the presence of the GHZ component 
being detected by state tomography. W states have been produced for studying 
disentanglement erasers [86j87f88] . 

Studies on larger spin system have mostly concentrated on cat states, which 
are analogous to GHZ states; notable results include 7-spin and 10-spin 
[10] states, as well as partial control of a 12-spin cat state [ID]- Cat states 
have also been explored in the context of entanglement assisted magnetic 
field sensing, where 10-spin |114] and 13-spin [115] states have been prepared 
with global control. The similarities between cat states and maximal quantum 
coherence have been explored in liquid crystal [132fl33|134fl35] and solid state 
systems. 



9 Conclusions 



Since the first tentative suggestions that NMR could be used to implement 
quantum computation [7] there has been quite extraordinary progress. In this 
final section I consider briefly how the field has stood up to early predictions, 
what has been achieved on the way, and where the future might lie. 

Since the very beginning of the field it has been understood that NMR can- 
not provide a practical technology for implementing genuinely useful quantum 
computations, partly because of the exponentially small signal size available as 
the size of the spin system is increased p^, but also because of the sheer diffi- 
culty of implementing quantum circuits in larger molecules [534] . The general 
belief was that systems with less than five qubits would be straightforward, but 
that larger systems would become increasingly difficult, with a limit between 
ten and twenty (just possibly thirty) qubits, and this has proved broadly cor- 
rect, with the largest spin-systems explored to date containing fourteen nuclei 

m- 

As expected [534] . NMR remained the leading technology for implementing 
quantum algorithms for around five to ten years. This reflects the extreme ex- 
perimental challenges involved in implementing quantum computations with. 
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for example, trapped ions. NMR remains by far the simplest technology for 
building small demonstration devices, and is so far the only technology used 
to implement Shor's quantum factoring algorithm. However the insurmount- 
able difficulties in building genuinely useful computers with liquid state NMR 
have led some authors to question whether there is any point in pursuing the 
field. If we cannot possibly reach our desired destination, what is the point 
in setting out? This question might seem reasonable, but its premise is too 
restrictive: sometimes one walks not to reach a destination, but to observe the 
scenery along the way, and the pursuit of NMR quantum computation has 
thrown up some surprising sights. 



Beyond the sheer interest of these ideas, one must also consider the potential 
role of NMR quantum computing in catalysing technology transfer between 
conventional NMR and quantum information processing. To date NMR has 
largely been the donor of ideas, with the use of composite pulses to suppress 
systematic errors being a particularly important example. Large scale quan- 
tum computers require single qubit rotations to be carried out with extremely 
high precision, and without the use of composite pulses it is hard to see how 
these could in fact be implemented experimentally. The quantum computing 
world has extended these ideas with the concept of arbitrarily accurate com- 
posite pulses [222], but it seems unlikely that these will find application in 
conventional NMR where ultra-high precision rotations are unnecessary. How- 
ever the near ubiquity of techniques to suppress pulse-length errors in NMR 
quantum computation should provide inspiration to more conventional work: 
the BBl composite pulse |213] is such a good replacement for a simple hard 
pulse that there is no reason why it should not be used completely routinely 
in most simple pulse-acquire experiments. 



Dynamical decoupling [255j is another area where the two fields may have 
much to offer one another, but it is perhaps not sensible to make too many 
firm and detailed predictions. Implementing quantum computing has enriched 
NMR by setting a series of challenging technical demands, which NMR has 
risen to magnificently. In return NMR has enriched quantum computation by 
proving that implementing complex unitary transformations can actually be 
quite straightforward, and by raising uncomfortable questions about the role of 
entanglement in mixed states. Most of the "easy" experiments have probably 
been performed by now, but there is still plenty to do, and the transfer of 
many of these experiments from liquid state NMR to more demanding fields, 
such as sohd state ENDOR where high purity states can be reached directly 
by cooling, is likely to keep us busy for some time to come. 
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